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Preface 

The  purpose  of  this  study  was  to  develop  a  computer  program 
which  would  provide  an  accurate  method  of  predicting  the 
stress/strain  behavior  of  advanced  composite  laminates.  The  ability  to 
predict  the  behavior  of  materials  at  large  strain  will  allow  the  use  of 
important  ply  lay-ups  such  as  ±45°  to  its  full  capacity. 

This  thesis  was  clearly  a  group  effort.  I  need  to  thank  my  thesis 
advisor.  Dr.  A.  N.  Palazotto,  for  his  patience  through  some  major 
problems.  I  am  also  deeply  indebted  to  Dr  R.  S.  Sandhu  of  the  Flight 
Dynamics  Laboratory  for  his  continued  assistance  over  the  past  year. 
Both  of  these  individuals  provided  me  with  an  insight  of  composites 
and  computer  programming  that  will  benefit  me  throughout  my 
career. 

The  people  of  the  structures  division: 

1.  Capt  J.  Daniels  and  Mr.  K.  Spitzer  for  insuring  the  specimens  were 
C-scanned,  cut,  and  ready  for  testing.  Capt  Daniels  also  helped  me  in 
understanding  the  computer  systems. 

2.  Mr  C.  Hitchcock,  and  the  instrumentation  group  who  attached  the 
strain  gages. 

3.  Mr.  D.  Cook,  Mr.  L.  Bates,  and  Mr  H.  Stalnaker,  of  the  Fatigue, 
Fracture,  and  Reliability  Group,  who  ran  the  testing. 

4.  Mr  Gene  Maddux  and  Don  Webb,  of  the  Photomechanics  Facility, 
also  provided  support. 
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Abstract 

A  geometric  nonlinear  technique  is  incorporated  in  a  current 
finite  element  program.  This  nonlinear  program  allows  material 
nonlinearity  for  calculating  the  stresses,  strains  and  failure  of 
composites.  The  improved  program  uses  an  updated  Lagrangian  to 
calculate  the  stresses  and  strains.  In  addition,  it  updates  the  fiber 
orientation  due  to  displacement  in  order  to  calculate  the  updated 
stiffness  matrix.  This  method  is  valuable  for  large  strain  values. 

The  analytical  data  were  compared  to  experimental  data 
obtained  from  Graphite  PolyEtherEther  Ketone  (Gr/PEEK)  [±45°]4s 

laminates.  To  obtain  data  two  geometries  were  used.  Digitized 
photographs  were  used  to  measure  the  angle  change  for  large  strains,  j  r  /  & 
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INVESTIGATION  OF  STRAIN  CHARACTERISTICS 
OF  GRAPHITE  POLYETHERETHER  KETONE 
USING  A  NONLINEAR  ANALYSIS 
AND  EXPERIMENTAL  METHODS 


L  Introduction 

Since  the  airplane  first  flew,  persistent  efforts  have  been  made 
to  make  it  lighter  and  stronger.  The  recent  development  of 
composites  has  added  to  this  impetus.  Graphite  based  composites  are 
stronger  and  lighter  than  aluminum.  Researchers  have  reduced  the 
weights  of  engines,  wings  and  several  other  aircraft  parts  through  the 
use  of  composites. 

In  the  past  several  years  new  composites  have  improved  aircraft 
performance  immeasurably.  The  Air  Force  uses  hundreds  of  different 
composites  in  aircraft  parts;  researchers  have  introduced  many  of 
these  materials  in  the  past  few  years.  Research  at  the  Flight  Dynamics 
Laboratory  of  Wright  Research  and  Development  Center  have  included 
many  of  these  composites. 

Graphite-epoxy  is  one  of  the  most  commonly  used  composites. 
The  materials  industry  has  recently  introduced  a  potential  substitute, 
graphite  polyetherether  ketone  (Gr/PEEK).  This  thermoplastic 
material  weighs  less,  has  a  greater  operating  temperature,  a  d  has 
better  fracture  toughnes,  characteristics  than  graphite-epoxy. 
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Composite  materials  are  the  macroscopic  blend  of  two  or  more 
different  materials.  Composites  include  cloth,  wood,  and  several  other 
materials  which  have  been  used  for  centuries.  This  research  is 
concerned  with  "advanced"  composites.  Advanced  composites  are 
materials  with  high  strength  fibers  blended  into  a  plastic  or  metal 
matrix. 

Manufacture  of  the  composite  can  taylor  differing  properties  as 
required.  When  the  composite  prepreg  is  manufactured,  the  fibers  are 
laid  such  that  the  direction  of  the  fibers  is  parallel.  Each  layer  or 
laminae  can  be  laid  with  a  different  orientation.  (Figure  1-1)  The 
matrix  is  then  allowed  to  cure  in  an  autoclave.  The  different 
temperatures  of  the  curing  process  can  produce  different  material 
properties.  The  data  sheet  from  the  production  of  the  composites 
used  in  this  thesis  is  provided  in  Appendix  A. 

"For  the  design  to  be  satisfactory  for  an  aircraft  it  is  essential  to 
determine  the  stress-strain  behavior  and  the  ultimate  strength  of  the 
laminates."  (21:104)  One  way  to  examine  the  stress-strain 
characteristics  is  through  computer  modeling,  using  finite  element 
models.  Another  way  is  by  experimental  observation.  The  best  way  is 
to  do  both  and  compare  the  results.  This  is  the  basis  of  this  thesis. 

A.  Problem. 

This  research  will  compare  the  stress-strain  analysis  of  several 
finite  elements  models  with  experimental  data.  Two  separate 
specimen  geometries  were  evaluated  numerically  .  (Figure  1-2)  The 
results  of  the  models  will  be  compared  to  experimental  data. 


I 
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Figure  1-1.  Usual  Fiber  Orientations 
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These  specimen  shapes  were  chosen  for  their  high  degree  of 
strain  provided.  Originally  a  0.4  inch  diameter  notched  model  was 
considered,  however,  to  obtain  high  strains  a  0.2  inch  diameter  notch 
was  used.  The  width  was  increased  to  1.125  inches  for  increased 
strain. 


1  1/8” 
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1  1/2" 
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I 

1  1/2' 
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Figure  1-2.  Specimens  Dimensions 
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Other  methods  of  obtaining  high  strain  results  are  discussed  in  a 
recent  article.  (11)  The  strain  field  mappings  similar  to  the  contours 
shown  in  Section  V. 

Earlier  versions  of  the  modified  program  were  used  on  several 
previous  theses.  (6,8,13)  This  thesis  will  update  the  program  PLSTR2 
to  provide  geometric  nonlinearity  along  with  the  material  nonlinearity 
shown  in  Capt  Daniels's  thesis.  (6).  The  geometric  nonlinearity  will 
consist  of  updating  the  Lagrangian  (shape  functions)  coordinates  and 
updating  the  fiber  orientation  (ply  angle)  during  each  increment  of 
loading.  A  ±45°  orientation  was  used  to  observe  the  fullest  use  of  these 
updates. 

The  experimentation  will  consist  of  strain  gage  measurements, 
and  extensometer  readings.  The  strain  gages  and  the  attaching 
adhesive  are  designed  for  up  to  10-20  percent  strain  per  the 
manufacturer.  Photographs  were  taken  during  the  loading  process  to 
determine  angle  changes  and  total  displacement  versus  loading. 

B.  Current  Knowledge. 

Three  recent  AFIT  theses  examined  the  characteristics  of 
Gr/PEEK.  (6,8,13)  Capt  Martin  studied  the  characteristics  at  room 
temperature  with,  and  without  concentric  holes  in  the  coupons.  Capt 
Fisher  studied  the  characteristics  at  high  temperature  with  and 
without  concentric  holes.  Capt  Daniels  studied  the  characteristics  at 
room  temperature  with  eccentric  and  concentric  holes.  Each  thesis 
used  different  methods  of  evaluating  the  stress-strain  characteristics. 

All  of  these  theses  looked  at  ±45°  laminates.  They  all  had 
difficulty  with  the  high  strains  obtained  by  these  laminates.  All  three 
theses  will  be  discussed  and  compared  with  the  present  research. 
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This  research  will  be  conducted  solely  on  (±45°)s  Gr/PEEK  16  ply 
laminates. 

These  theses  also  contained  finite  element  models.  Capt  Daniels 
used  the  program  called  PLSTREN  (6:6-4)  which  incorporated 
nonlinear  material  properties.  This  research  will  take  that  same 
program  and  account  for  geometric  nonlinearity  as  well.  These 
changes  are  discussed  fully  in  Section  III. 

This  type  of  analysis  is  necessary  for  large  strain.  Mr.  Nabil  Y. 
Ghantos  wrote  a  program  for  "Geometric  Nonlinear  Analysis  of  Plane 
Frame  Structures".  His  thesis  shows  the  enhanced  predictions  of 
using  geometric  nonlinearity  in  plane  structures.  (9,10-19)  His  thesis 
shows  that  if  structures  experience  large  stresses  the  shape  of  the 
structure  will  change.  Therefore,  during  the  loading  process  the 
change  of  new  shape  of  the  structure  must  be  updated.  This  thesis 
uses  that  same  concept  only  with  a  composite  coupon  instead  of  a 
structure. 

Large  strains  have  always  been  a  difficult  process  to  predict 
since  the  process  is  highly  nonlinear.  Linear  equations  are  not 
accurate  at  high  strains,  and  so  the  nonlinear  analysis  is  necessary.  All 
three  AFIT  theses  done  on  this  material  looked  at  ±45°  lay-ups.  Most 
of  the  material  properties  which  Capt  Martin  (13:116-118)  derived 
will  be  used  in  this  thesis.  However,  he  decided  to  select  a  cut  off 
point  of  5%  strain  to  separate  interlaminar  and  intralaminar  ahear 
strains.  By  using  high  elongation  gages  and  high  elongation  adhesive, 
much  higher  strains  than  either  Martin  or  Fisher  were  observed. 
Therefore,  shearing  stress  and  strain  will  be  more  accurate  in  these 


tests. 


Several  articles  have  been  written  on  the  subject  of  notched 
laminates  (20).  Notched  specimens  were  also  looked  at  in  this  study 
and  will  be  compared  with  previous  investigations. 


IL  Theory 


The  applicable  theories  used  in  this  thesis  include  the 
constitutive  equations  of  composite  materials,  failure  theory,  and 
nonlinear  finite  element  theory,.  The  constitutive  laws  are  utilized  in 
the  program  and  in  the  analysis  of  the  data  obtained  experimentally. 

This  research  will  be  looking  at  ±45°  ply  lay-ups.  This  will 
create  large  displacements  and  thus  strains  and  the  fiber  orientation 
will  change  as  the  specimen  strains.  It  is  assumed  that  plasticity  will 
not  occur,  but  stress-  strain  relations  can  be  considered  nonlinear. 
Figures  2- la,  b  and  c  show  how  the  lay-up  starts  and  is  updated  as 
stresses  are  applied.  These  figure  show  the  directions  of  fiber 
orientation.  Figure  2-lc  show  how  the  angle  changes  as  the  forces  are 
applied.  The  stresses  in  the  X  direction  are  derived  from  the  applied 
forces.  The  stresses  in  the  Y  direction  are  zero.  The  stresses  in  the 
fiber  directions  change. 

The  basic  constitutive  equations  include  Hooke's  Law  equation  2- 
1.  (18:1881 

°i  =  Cijej  i  =  1.-6  j  =  1..6  (2-1) 


where: 

°i  =  stresses 

Cy  =  stiffness  matrix 

ej  =  strains 
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Figure  2- la.  Stress  Directions  in  -45°  Specimen 


Figure  2- lb.  Stress  Directions  in  +45°  Specimen 


Figure  2-lc.  Stress  Directions  After  Deformation 

For  our  purposes,  only  plane  stress  and  plane  strain  will  be 
looked  at.  In  other  words,  three  dimensional  effects  will  not  be 
considered.  a3,T13.T23=0:  and  £3 ,^13,^23=0 

Equation  (2-2)  shows  the  transformation  matrix  used  to  solve  for  the 
stresses. 


dox 

-2mn 

dax 

<  da2 

>  = 

m^  2mn 

< 

day  ” 

.dx12. 

-  run  -mn  m2-n2  - 

k  dxxy  > 

(2-2) 


Where: 

d^x  =  change  of  stress  in  direction  of  specimen, 
doj  =  change  of  stress  in  fiber  direction 
m  =  cos  0 
n  =  sin  9 
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0  =  fiber  orientation  (updated  throughout  program) 


A.  Linear  Finite  Element  Theory. 

Finite  element  analysis  utilizes  the  equations  described  above.  It 
assumes  that  each  element  has  the  same  properties  throughout  the 
element.  At  points  where  there  is  a  stress  concentration,  a  closely 
refined  model  is  needed.  In  our  study  as  many  elements  as  feasible  to 
get  the  best  results  will  be  considered.  Figure  2-3a  shows  the  original 
model.  Figure  2-3b  shows  the  quarter  model.  Symmetry  was  taken 
advantage  of  in  the  model. 


Figure  2-3a.  Original  Half  Model 
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Figure  2-3b.  Quarter  Model 

The  program  PLSTR2  uses  quadrilateral  elements  consisting  of 
constant  strain  triangular  elements.  This  concept  is  described  in 
detail  in  several  books  (2:98-101).  The  use  of  the  concept  in 
PLSTREN  (linear)  and  PLSTR2  (nonlinear  material)  is  described  in 
detail  by  Capt  Daniels.  (6:2-15  -  2-19)  The  constant  strain  triangle  is 
shown  in  figure  2-9. 

All  element  in  this  program  are  "Lagrange  elements"  (2:97-99). 
Therefore,  the  elements  shape  functions  are  derived  using  "Lagrange's 
Interpolation  Formula".  This  formula  is  reduced  to  (2-4)  for  plane 
stress/plane  strain  problems. 

n 

u  =  XNjUj  or  u=N1ui+N2u2+ . +  Nnun  (2-4) 

i=l 

n 

v  =  XNiV;  or  v  =  NiVi  +N2v2  + . +  Nnvn 

i=l 


where:  N=  shape  function 
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u=  X  direction  component  of  displacement 
v=  Y  direction  component  of  displacement 
n=  4  for  a  four  noded  element  (3  in  this  program) 

This  program  will  divide  all  elements  into  triangles  as  discussed 
further  on  in  this  section. 

The  models  used  in  this  thesis  use  both  triangular  elements  and 

condensed  rectangular  elements.  The  elements  must  be  kept  as  close 

to  square  as  possible  close  to  either  a  material  or  geometric 

discontinuity.  (Figure  2-4).  This  program  updates  the  coordinates  and 

thus  the  element’s  shape  is  updated  as  shown.  Triangular  elements 

are  only  used  at  uniform  stress/strain  locations  (Figure  2-5). 

Original  Shape  Updated  Shape 


Figure  2-4.  Degenerated  Rectangular  Elements 


Figure  2-5.  Triangular  Elements 

Elements  are  numbered  counterclockwise.  The  centroid  is  then 
calculated  and  all  stresses  and  strains  are  calculated  for  that  point. 
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These  relationships  are  used  in  both  linear  and  nonlinear 
programs.  The  linear  program  is  only  used  to  compare  with  the 
nonlinear  program  and  show  the  integrity  of  our  models. 

B.  Material  Nonlinear  Constitutive  Relations. 

The  material  properties  of  Gr/PEEK  have  been  shown  by  Capt 
Martin  to  change  as  the  material  is  deformed.  Material  properties 
were  obtained  experimentally  for  tension,  compression,  and  shear.  He 
averaged  the  data  by  selecting  suitable  sets  of  strain  values  for  all  the 
stress  strain  curves  required  in  this  computer  code. 


Figure  2-6.  Examples  of  Cubic  Splines  Used  in  Program  PLSTR2 
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As  each  increment  of  load  is  applied,  new  material  properties 
are  calculated  using  a  predictive,  corrective,  and  iterative  technique. 
This  process  is  described  in  detail  by  Capt  Fisher.  (8:2-15  -  2-29) 

After  enough  load  is  applied  to  extend  past  the  end  of  the  curve, 
the  program  believes  the  material  has  failed.  For  our  models  failure 
occurred  at  approximately  6  times  the  level  of  strain  experienced  by 
Capt  Martin.  Therefore,  experimental  data  from  these  experiments 
were  used  for  the  shear  stress\strain  curve  with  angle  updates.  (Figure 
4-6).  The  curves  for  glass  epoxy  were  also  extended.  (Figure  2-7) 


Strain  in  fiber  direction  (in/in) 

Figure  2-7  Extended  Curve  of  Capt  Martin's  Data 


C.  Geometric  Nonlinear  Constitutive  Relations. 
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The  geometric  nonlinearity  provided  by  this  program  was  done 
in  two  parts.  The  coordinates  were  updated  after  every  increment  of 
load  was  applied.  Each  time  the  strain  increments  were  calculated 
the  fiber  orientation  of  each  ply  was  updated.  A  description  of  the 
updated  Lagrangian  coordinate  system  is  given  by  Owen  and  Hinton 
(15:382-388) 

To  update  the  coordinates,  one  first  must  calculate  the  Q  matrix. 
By  multiplying  the  Q  matrix  (equation  2-4)  by  the  transformation 
matrix(equation  2-1)  and  it's  inverse  one  can  obtain  the  Q 
matrix(equation  2-4). 


r  > 

dox 

Q11  Q 12  0 

r  dc. 

•< 

da2 

>  = 

Ql2  Q22  0 

< 

dC2 

(2-4) 

-  0  0  Q66  _ 

-  de3  . 

Where: 

Q11  =  - 

Ql2= 

- 

Ejv12 

1-v12v21 

'-Wll 

Q22  =  7  ~ —  Q66=Gi2 

1-Vi2v2i 


r  > 

dax 

Q11  Qi  2  Qi  6 

r  ^ 

dEx 

*  doy 

►  = 

Qi  2  Q22  Q26 

< 

dEy  ^ 

dxvv 

v  xy  j 

-  Oi 

.  dTxy  , 

Where:  Qi  1  =  Qi  1m4  +2(Qi2+2Q66)n2m2  +Q22n4 

Ql2=  (Qn+Q22-4Q66)n2m2  +Q22(n4+m4) 
Q22  =  Qnn+2(Qi2+2Q66)n2m2+Q22m4 


I 
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Qi6  =  (QirQi2-Q66)nm3  +  (Qi2-Q22+2Q66)n3m 

Q26  =  (QirQl2-Q66)n3m  +  (Ql2*Q22+2Q66)nm3 

Q66=(Qii+  Q22-2Qi2*2Q66)n2m2  +  Q66(n4+m4) 

m  =  cos0 
n  =  sin9 
0  =  fiber  orientation 

Once  the  Q  matrix  has  been  calculated,  one  can  determine  the 
incremental  strains  from  the  incremental  stresses.  From  the 
incremental  strains  one  can  now  calculate  the  displacements. (equation 
2-8)  These  equations  were  shown  by  Palazotto  (16). 


[Keq]  =  aLbFIQHBH,  (2-6) 

!  =  1 

n 

[dN]  =  X  lQ]k  ldulk  lk  (2-7) 

p=i 

[dul  =  (Keq]_1[dN]  (2-8) 

where:  k  =  number  of  plies 

=  thickness  of  kth  ply 

After  the  displacements  are  calculated,  the  new  element  is 
defined.  One  now  has  a  new  mesh  from  which  to  calculate  stresses 
and  strains  (Figure  2-8). 
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Figure  2-8  a.  Old  Mesh 


Figure  2-8  b.  New  Mesh 

Calculation  of  the  new  ply  orientation  is  done  prior  to  calculating 
the  Q  matrixfequation  2-9).  The  accuracy  of  this  equation  is  shown  in 


figure  5-7  of  the  results. 

d0  =  dYxyUS0/71) 


Where: 


d0  =  change  in  fiber  orientation(degrees) 
Yxy  =  shear  strain 


(2-9) 


n  Program  Description. 


'The  laminate  constitutive  relations  are  initially  assumed  to  be 
linear  and  the  element  stiffness  matrices  are  computed"  (19:182) 

This  stiffness  matrices  are  made  into  a  global  stiffness  matrix  and  then 
condensed  into  a  reduced  stiffness  matrix.  The  resulting  equations 
are  solved  using  the  applied  incremental  loads.  These  equations  are 
now  used  to  solve  for  the  displacements. 

The  calculated  displacements  from  an  incremental  point  of  view 
can  be  written  as  vectors  using  the  updated  Lagrangian  approach. 

The  updated  Lagrangian  is  acceptable  from  the  small  increments  of 
displacement. 

du=[X]{ai  a2  a3}T  and  dv  =(X]{a4  as  aeP  (2-10) 

where:  [Xl=[l  xy]  from  figure  2-7 

x  and  y  are  the  coordinate  differential  functions 

These  displacements  are  converted  into  strain.  The  higher 
order  terms  of  the  total  Lagrangian  are  ignored  due  to  small 
increments. 


ddu 


ddv 


ddu  .  ddv 


ip  uuu  i r*  vuv  iai  vuu  i  w v-*  i 

d8x=^r  d£y*air  dyxy=F^ 


(2-11) 


These  strains  are  now  used  to  calculate  the  strain-displacement 
matrix  [B] . 


r 

d£x 


V. 


d£y 

dYxy 


M  BIS 


(  dui") 
dvi 

du2 
dv2 
du3 
^  dv3 


(2-12) 
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where: 

-  1  0  0  0  0  0  - 

[B]=  0  0  0  1  0  0  [A]1 

-  1  0  1  0  0  0  - 

"  1  xl  yl  “ 

[A]=  1  x2  y2  previous  coordinate  before  the 

-  1  x3  y3  - 

increment  of  displacement. 

[A]  is  calculated  from  the  constant  area  triangle  shown  in  figure 
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Figure  2-9.  Constant  Strain  Triangle 

The  stiffness  matrix  Ik]  of  an  element  can  now  be  calculated 
from  equation  2-13. 

[k]e=JA  [BF  IQ]|B|tdA=[BF  [Q]|B]tA  (2-13) 

where:  A  is  previous  element  area 
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At  previous  load 


Q  = 


Qii  Qi  2  Qi  6 
Ql  2  Q22  Q26 
“  Ql  6  Q26  066  - 


Once  the  new  strains  are  calculated,  and  the  position  of  the 
nodes  are  determined,  with  these  new  positions  it  is  possible  to 
evaluate  updated  element  areas  and  a  check  to  see  if  the  maximum 
strain  has  been  reached  for  ply  failure.  It  is  then  possible  to  calculate 
the  stresses  by  adding  the  new  increment  of  stress  to  previous  total 
stress.  A  more  complete  discussion  of  these  equations  is  provided  in 
section  III. 

The  failure  criterion  used  on  PLSTREN  and  PLSTR2  was 
discussed  in  many  references  and  is  explained  fully  in  Technical 
Report  AFFDL-TR-73- 137  (23). 

For  this  program,  total  strain  energy  determines  the  failure  of 
laminates. 

*1*1  +  K^2  +  KgWg  >  1  (2-14) 


Where: 

Wi  =Je,  Ojdej  1/IS  =leju  ot dej  (2-15) 

The  program  considers  the  laminate  as  having  failed  if: 
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KlWl/(KlWl  +k2w2  +  k6w6) 


>0.1  for  fiber  failure 

(2-16) 

<0.1  for  matrix  failure 

Two  different  types  of  unloading  have  been  considered  for  this 
program.  However,  gradual  unloading  is  most  accurate  for 
multidirectional  laminates.  For  this  program  this  is  accomplished  by 
setting  NOPSHN  equal  to  2.  For  immediate  unloading,  set  NOPSHN 
equal  to  1.  Upon  failure  the  modulus  of  elasticity  for  the  failed 
laminate  is  made  negative.  This  unloads  the  laminate  as  the  load 
increases  and  has  been  shown  to  be  quite  accurate  with  experimental 
data.  (20:168) 

This  failure  criteria  is  compared  to  numerous  other  criteria  by 
Sandhu.  (22)  This  method  shows  excellent  results  when  compared  to 
experimental  data. 
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UL  Analysis 

This  chapter  includes  a  description  of  the  geometry  of  the 
specimen,  how  the  models  were  developed,  and  a  convergence  study 
to  determine  the  integrity  of  the  models. 

This  thesis  compares  the  linear  finite  element  model  with  the 
nonlinear  models.  The  program  diagram  has  not  changed  from  Capt 
Daniels's  thesis.  We  follow  the  subroutines  on  the  nonlinear  part  of 
the  program.  The  changes  in  the  program  were  in  the  STON,  and 
OUTPUT  subroutines. 

This  thesis  uses  a  conventional  method  of  Finite  Element 
Analysis  (FEA).  Betts  discusses  several  different  methods  of  FEA  and 
commercially  available  FEA  programs.  As  Betts  states,  the  re-mesh  of 
the  model  may  create  a  potential  for  distorting  element  shapes(l:60). 
This  could  decrease  accuracy.  To  show  the  integrity  of  our  model  we 
will  compare  with  experimental  results.  The  convergence  study  also 
proves  the  integrity  of  the  models. 

A.  Finite  Element  Modeling. 

This  thesis  modeled  the  specimens  shown  in  Figure  1.2.  The 
modeling  consisted  of  highly  refined  elements  around  the 
discontinuities.  This  included  both  material  and  geometric 
discontinuities.  Elements  in  between  were  made  as  geometrically 
convenient  since  stress  and  strain  were  reasonably  constant.  The  use 
of  triangles  was  a  mistake.  The  program  works  better  with  four  sided 
polygons  as  shown  by  Sandhu(20).  For  this  reason  a  smaller  model 
was  used  for  some  of  our  calculations.  As  shown  in  the  section  V,  the 
smaller  model  provided  better  failure  criteria.  (Figure  5-4c) 
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The  material  properties  for  glass  epoxy  the  linear  program  were 
taken  from  Cron's  experimental  results.  (3)  Cron  ran  tests  on  both 
Gr/epoxy  and  glass  epoxy.  His  data  had  a  constant  value  for  the 
material  properties.  These  properties  of  glass  epoxy  were  also  used  by 
Capt  Martin. 

Using  Capt  Martin's  model,  it  was  observed  that  the  differences 
in  using  the  quarter  model  were  small. (Figure  3-1)  This  graph  shows 
the  differences  between  strains  from  the  linear  program  across  the 
line  of  symmetry  (Figure  3-2).  Since  the  differences  were  less  than  1 
percent,  we  considered  the  model  viable. 


3-2 


Distance  From  the  Center  of  the 
Discontinuity  (in) 


Figure  3-1.  Half  and  Quarter  Model  Differences 


Therefore,  we  assume  that  the  differences  are  minimal  using  our 
models  for  determining  stresses  along  the  line  of  symmetry.  (Figure  3- 
2)  All  other  stresses  and  strains  were  nearly  the  same.  Figure  3-3 
shows  the  model  used  herein.  This  model  has  629  nodes  and  1284 
elements  which  is  about  half  the  nodes  and  elements  of  the  half 
model. 


3-3 


Difference  along  line  of  symmetry  0  6-  from  center 


Figure  3-2.  Size  of  half  and  Quarter  Models 


Figure  3-3.  Quarter  Notched  Finite  Element  Model 

The  area  around  the  discontinuity  was  modeled  using  a 
program  to  divide  the  radius  of  the  notch  by  a  specified  number.  This 
was  the  size  of  sides  of  the  elements  nearest  the  notch.  From  there, 
the  elements  were  expanded  very  slowly  to  create  a  highly  refined 
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model.  The  minimum  size  of  the  elements  around  the  notch  is 

0.0046"  by  0.0046"  and  were  square.  The  maximum  size  is  0.5"  by 

0.1"  and  has  an  aspect  ratio  of  five. 

Side  View  of  Specimen 


Tab  Taper 


Glass  Epoxy 


16  plys 
Gr/PEEK 


Figure  3-4.  Side  View  of  Modeled  Region 

The  actual  lay-up  of  the  (±45°]s  laminate  has  16  plies.  This 
program  assumes  only  two  plies.  Each  ply  is  8  times  as  thick  as  one 
actual  ply.  In  other  words,  instead  of  modeling  16  .00525"  plies,  eight 
.042"  plies  were  modeled. 


B.  Modeling  of  the  Glass  Epoxy  Fibers. 

As  shown  in  figure  3-4,  the  taper  of  the  tabs  was  modeled.  The 
material  properties  from  Capt  Martin's  thesis  were  used  for  the  lower 
strains.  However,  the  cubic  spline  data  was  extrapolated  up  to  40 
percent  strain.  Although  Capt  Daniels's  thesis  modeled  the  entire  tab, 
it  was  not  necessary  for  this  thesis,  since  the  fixture  is  assumed  to 
remain  undeformed  and  provide  load  according  to  the  load  cell. 

C.  Convergence  Studv. 


The  convergence  study  was  done  in  the  same  way  as  Capt 
Daniels,  and  Capt  Fisher,  These  were  only  quarter  models  of  the 
specimen,  so  we  needed  to  do  a  study.  This  is  done  by  running  the 
model  as  an  isotropic  material  using  the  linear  finite  element  program. 
From  Peterson  (17:21-24)  and  Griffel  (10),  equation  (3-1)  is  applied. 
°max 

Kt  =  -  (3-1) 

Onom 

Kj.  =  Stress  Concentration  Factor 
a  =  load  divided  by  area 

Peterson  and  Griffel  both  give  Kt  =3.065  for  a  circular  notch. 

Using  the  material  properties  from  Capt  Martin,  our  notched  model 
gives  Kt  =  2.915.  The  notched  model  decided  on  was  within  4.89% 

which  was  adequate.  This  is  not  surprizing  since  Capt  Daniels's  model 
with  a  hole  had  a  difference  of  4.64%  and  Capt  Fisher  4.93%. 

For  the  case  without  any  discontinuity  there  is  no  stress 
concentration  or  Kt  =  1.0.  For  the  small  model  (Figure  3-5)  the  stress 

concentration  varies  up  to  1.029  or  2.9%  error.  The  large  model 
difference  is  small  enough  to  not  show  up  on  the  computer  print  out. 

Part  of  the  reason  for  the  error  is  that  a  quarter  model  was  used. 
A  half  model  would  have  a  better  ratio. 
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Figure  3-5.  Small  Quarter  Model  No  Discontinuity 
D.  Boundary  Conditions 

The  boundary  conditions  are  shown  in  Figure  3-6.  Nodes  along 
the  y-axis  were  fixed  in  the  x  direction.  Nodes  along  the  x-axis  were 
fixed  in  the  y  direction.  Nodes  along  the  far  side  of  the  x-axis  were 
displaced  incrementally.  All  other  nodes  were  allowed  to  move  freely 
in  either  direction. 
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No  Discontinuity 


Notched  Specimen 


A 

Figure  3-6.  Boundary  Conditions 
This  allows  for  2  degrees  of  freedom  for  the  unrestrained  nodes. 
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IV.  Experimentation 


The  purpose  of  the  experimentation  in  this  thesis  is  to  compare 
the  results  with  the  analytical  models.  Experimentation  was 
conducted  at  the  facilities  of  the  Structures  Division,  Flight  Dynamics 
Laboratory  (FDL),  Wright-Patterson  AFB,  Ohio.  Set  up  is  shown  in 
Figure  4- 1 .  The  voltmeter  was  used  to  show  the  load  at  the  time  the 
photograph  was  taken.  For  Figure  4- lb,  take  the  voltage  minus  the 
voltage  with  no  load  and  divide  by  two.  This  gives  us  510  lbs  applied 
at  the  time  the  photograph  was  taken.  During  the  testing  this  number 
was  checked  with  the  value  of  the  load  cell  and  found  to  be  almost 
identical. 


Figure  4- la.  Photograph  of  Setup 
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Figure  4- lb.  Close-up  Photograph  of  Setup 

The  voltmeter  is  used  to  show  the  load  applied.  To  calculate  the 
load,  subtract  the  initial  voltage,  divide  by  two,  and  multiply  1000.  For 
figure  4- lb,  the  load  is  510  pounds. 

The  two  models  shown  in  figure  4-2  were  tested  using  seven 
specimens  of  each  type.  High  elongation  gages  were  used  for  all 
experimentation.  Stacked  rosettes  were  not  available  for  high  strains. 


therefore,  rosettes  are  per  figure  4-3.  A  complete  list  of  materials 

used  is  provided  in  Appendix  B 

Gages  located  in  center  of  specimen  and  1“  from  tabs 


Inked 

Cross 


Back  to 

Back 

Rosettes 


Rosette 


Radius  =  0. 1  inch 


Inked 

Cross 


Back  to 
Back 
Rosettes 


Rosette 


Figure  4-2.  Gaged  and  Crossed  Specimens 


Strain 


Figure  4-3.  Electrix  Industries  High  Elongation  Rosettes 
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The  Gr/PEEK  panels  were  manufactured  by  the  Fiberite 
Corporation,  a  subsidiary  of  Imperial  Chemical  Industries  of  Great 
Britain.  The  production  information  is  described  in  the  data  sheets 
(7).  The  purchase  order  is  provided  in  Appendix  B.  These  panels 
were  C-scanned  to  check  for  preexisting  flaws.  None  were  found. 
(The  C-scan  is  provided  in  Appendix  B)  Only  (±45°)s  16  plies  were 
used. 

Table  4-1  Specimen  Dimensions 


Laminate  Type 

Ply  Lay-up 

Size 

Ply  Thickness 

No  discontinuity 

(±45)s 

10"  X  1" 

16  plies 

Notched  specimen 

(.1"  radius) 

(±45)s 

10"  X  1.2" 

16  plies 

After  C-scanning  the  panels  were  cut  into  different  sizes  shown 
in  table  4-1.  Tabs  were  bonded  to  each  specimen  to  allow  the  Instron 
machine  to  grip  the  specimen.  The  tabs  were  1/16  inch  thick  G-10 
glass  epoxy  (0/90  woven)  and  were  bonded  to  the  specimens  using  an 
adhesive  which  required  curing  in  an  autoclave. 

High  elongation  strain  rosettes  (PAHE  -  03  -  125RB  -350  LEN 
and  PAHE  -  03-062RB  -  350  LEN)  were  used.  The  gages  were  tested 
for  proper  resistances.  Several  gages  were  not  within  the 
specifications  and  needed  to  be  replaced.  These  gages  were  attached 
with  Micro  Measurements  M-Bond  610  for  high  elongation.  These  are 
the  highest  gages  and  adhesive  found  for  this  type  of  experiment.  This 
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provided  us  with  higher  strain  data  than  Fisher.  The  gages  did  not  fail 
until  almost  30%  strain  was  reached. 

Since,  the  material  properties  from  Martin’s  experiments  were 
used,  no  compression  testing  was  done. 

Camera  work  was  provided  by  the  Technical  Photographic 
Division  (4950TESTW/RMP).  Pictures  were  taken  of  an  "X"  placed  on 
the  side  opposite  the  strain  gage.  This  showed  the  change  of  the  45° 
ply  angle  as  the  load  was  applied.  The  pictures  were  blown  up  and  the 
angle  digitized.  Figure  4-7  shows  the  change  in  angle  as  the  load  was 
applied  for  the  specimens  without  discontinuities.  This  was  the 
average  of  two  groups  of  photographs  taken,  one  per  specimen.  Figure 
4-8  shows  the  change  in  angle  for  the  notched  specimen.  Only  one 
set  of  photographs  were  taken. 

A.  Test  Apparatus. 

The  Instron  machine  used  to  run  these  tests  had  a  maximum 
load  capacity  of  20  kips.  It  was  run  at  0.05  inches  per  minute.  Strain 
gage  readings  were  recorded  once  every  four  seconds.  This  provided 
plenty  of  data,  since  each  specimen  with  no  discontinuity  took  45-50 
minutes  to  fail  and  each  notched  specimen  took  20-25  minutes  to  fail. 

A  data  file  was  created  on  the  Fight  Dynamics  Laboratory  (FDL) 
VAX  which  contained  all  the  gage  data,  the  load  data  and  the  time. 
From  this  a  small  program  was  created  to  divide  the  data  into  smaller 
sections  and  into  a  format  readable  by  Cricket  Graph,  the  graphics 
program  used  to  make  the  graphs  for  this  thesis.  By  sampling  only  50- 
60  points  out  of  2000-4000  we  only  missed  the  small  changes 
provided  by  the  fibers  failing  at  the  beginning  of  the  experiment. 
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Photographs  were  taken  at  approximately  2  feet  from  the 
specimen.  The  photographs  were  only  taken  of  three  different 
specimens.  Preliminary  photographs  were  taken  on  several  of  the 
other  specimens. 

B.  Preliminary  Results. 

Material  properties  from  Capt  Martin  and  Fisher’s  experiments 
were  valid  only  to  approximately  5-10  percent  strain.  The  present 
experiments  provided  better  data  at  the  higher  strains.  By  using  Capt 
Martin’s  finite  element  model  data,  the  program  predicts  failure  at 
slightly  more  than  .25  inches  for  a  specimen  with  no  discontinuity. 
Experimentally,  however,  failure  was  found  at  around  2.1  inches.  This 
is  due  primarily  to  the  data  being  good  only  to  about  5  percent  strain 
for  the  gages  used  by  Martin.  (Figure  4-4)  Since  the  angles  are 
changing,  the  shear  stress  and  strain  must  be  updated.  The  shear  data 
is  calculated  from  equations  4-1  and  4-4.  The  tabular  data  is  provided 
in  Table  4-2. 
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30000 


Shear  Strain  (micro  in/in) 


Figure  4-4.  Experimental  Shear  Data  Comparisons 

The  new  values  for  shearing  strain  were  calculated  using  equation  4-1 
and  4-2. 

Yl2  =  2(ex_ey)Sin0C°Se  (4-1) 


where:  Sc  =  axial  strain 

*y  =  transverse  strain 

0  =  fiber  orientation  (updated  throughout  program) 
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This  equation  comes  from  using  the  transformation  matrix 
shown  in  equation  4-2.  These  basic  equations  come  from  several 


sources.  (5:Section  2.1.2;  1-11)  (12:48-51) 
-  ^ 


*2 


Yl2 

v ~2~ J 


I-  m2  n2 
m2 


•2mn 

2mn 


>-  nm  -mn  m^-n^  _i 


Ex 
®y 

Yxy 

v  2  ^ 


(4-2) 


where:  m  =  cos  0 

n  =  sin  0 

0  =  fiber  orientation  (updated  throughout  program) 


y12  =  2(ex  -  Ey)sin0cos0  +  y^cos2©  -sin20) 

The  shear  strain  (y  j  is  small  because  the  rosette  is  located  at  the  line 
of  symmetry.  Therefore,  equation  4-1  was  used,  is  given  from  the  C 

leg  of  the  rosette.  £y  is  given  from  the  A  leg  of  the  rosette. 

The  same  transformation  matrix  can  be  used  to  calculate  shear 
stress  also. 


°1 

n2 

-2mn 

<  C2 

>  = 

n2 

m2 

2mn 

< 

°y  ' 

.T12, 

-  nm 

-mn 

m2-n2  - 

-  ^xy  - 

(4-3) 


where:  m  =  cos  0 

n  =  sin  0 


Assuming  txy  is  small,  equation  4-4  results.  This  thesis  does  not 
look  at  changes  to  the  width  or  thickness 
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cos©  sin© 


(4-4) 


*12  = 


P  =  load  at  end  of  specimen  (lbs) 
b  =  width  of  specimen  (1.0") 
d  =  thickness  of  specimen  (.084") 

0  =  fiber  orientation  (degrees) 

The  angle  was  calculated  using  the  cross  on  the  specimen. 
(Figure  4-5)  The  cross  was  placed  one  inch  from  the  tab,  the  same 
location  as  the  strain  gage,  only  on  the  other  end. 


Figure  4-5.  Angle  Calculation  from  Cross  on  Specimen  Displacement 
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Figure  4-6.  Photograph  of  Cress  on  Specimen 


Table  4-2  Shear  Uata 


Shear  Strain 
0.000 
0.005 
0.018 
0.034 
0.051 
0.070 
0.108 
0.129 


Shear  Stress 
0.0 
2986.0 
7743.0 
9552.0 
10420.0 
11205.0 
12730.0 
13517.0 
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0.170 

15146.0 

0.209 

16831.0 

0.248 

18530.0 

0.283 

20339.0 

0.318 

21724.0 

0.351 

23186.0 

0.378 

24690.0 

The  assumptions  used  by  Fisher  "...fibers  retain  their  ±45 
orientation."  (8,4-7)  This  was  admittedly  a  bad  assumption,  but,  at  low 
strain  rates  the  error  introduced  was  small.  Obviously,  this  is  not  a 
good  assumption  for  our  tests.  From  the  photographs,  one  can  see 
that  the  angle  change  is  significant.  (Figure  4-7)  The  angle  change  for 
the  notched  specimen  follows  a  similar  pattern.  (Figure  4-8)  The 
increased  load  prior  to  fiber  change  is  due  to  the  wider  specimen. 

Angle  changes  for  the  notched  specimen  showed  a  similar 
pattern.  It  can  be  seen  that  the  angles  show  little  change  before  two 
thousand  pounds  are  applied. 
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One  problem  encountered  was  the  strain  gages.  At  first  the 
strain  gages  saturated  at  10  percent  strain.  The  gages  saturated 
because  of  the  amount  of  current  applied.  This  problem  was  solved  by 
reducing  the  "gain  factor"  to  20.  The  gain  factor  is  described  in  detail 
by  a  gage  manufacturer  (14).  By  reducing  the  gain  factor  one  reduces 
the  amount  of  current  flowing  through  the  gage.  When  the  current 
gets  to  a  certain  level  the  gage  saturates.  This  reduced  the  accuracy  of 
the  readings  but  allowed  us  to  get  reading  beyond  10  percent  strain. 
(Figure  4-2)  The  gages  also  peeled  off  the  specimen  at  high  strains. 
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Figure  4-8.  Angle  Changes  in  Model  with  Notch 


This  problem  was  solved  by  running  the  experiments  as  soon  as 
possible  after  the  gages  had  been  mounted.  In  addition,  the  leads 
were  soldered  directly  to  the  specimen.  The  gages  pulled  off 
primarily  because  of  the  differences  in  axial  and  transverse  strain  at 
the  location  of  the  gages. 

The  tab  material  debonded  (Figure  4-9).  This  is  caused  by  the 
high  stress  concentrations  at  the  corners.  This  was  also  experienced 
by  the  other  three  theses.  During  the  high  levels  of  strain  it  is  more 
evident.  This  debonding  begins  at  approximately  5  percent  strain. 

The  only  way  to  alleviate  this  problem  would  have  been  to  increase  the 
length  of  the  specimen.  Since  the  plates  were  already  manufactured 
and  cut,  this  was  not  a  possibility. 
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Figure  4-9.  Sketch  of  Debonding 

A  stress/strain  graph  is  shown  in  Figure  4-10.  This  graph  shows 
the  high  rates  from  the  gages  without  failure.  This  data  is  taken  from 
the  average  of  the  two  central  gages  for  the  last  two  specimens. 
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Figure  4-10.  Stress  vs  Strain:  No  Discontinuity 
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V.  Results  and  Discussion 


This  chapter  will  compare  both  the  experimental  data  and  the 
analytical  data.  This  is  done  in  two  steps;  comparison  of 
experimentation  to  analytical  data  and  comparison  of  linear  to 
nonlinear  programs. 

A  Load.  Displacement  Comparison 

First  of  all.  plotting  the  standard  load  vs  displacement  curve 
shows  that  our  model  is  reasonably  accurate.  Figure  5- 1  shows  the 
load  versus  displacement  curve.  This  is  the  average  of  seven  tests 
done  on  specimens  with  no  discontinuity.  As  shown  is  section  IV,  the 
angle  change  is  only  necessary  at  high  displacements.  The 
experimental  data  was  taken  directly  off  the  extensometer  and  the 
load  cell.  (Figure  5-1)  For  the  analytical  models  we  used  an  average  of 
ox  along  the  edge  in  equation  5-1. 

P  =  ox  b  d  (5-1) 

P  =  load  at  end  of  specimen  (lbs) 
b  =  width  of  specimen  (1.0") 
d  =  thickness  of  specimen  (0.084") 

The  displacements  and  loads  were  taken  from  the  models 
shown  in  Figure  5-2.  Because  this  is  a  quarter  model  we  must  multiply 
displacements  at  the  edge  by  two.  The  nonlinear  programs  are  run  at 
displacement  increments  of  .01".  For  the  linear  program  the 
displacement  is  specified  at  the  start  of  the  program  and  run  at 
different  displacements. 
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Figure  5-1.  Experimental  Displacement 
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Figure  5-2.  Analytical  Displacements  and  Loads 
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Figure  5-3.  Load  vs  Displacement  Curve  for  No  Discontinuity. 

From  Figure  5-3,  the  updated  model  follows  the  experimental 
data  most  closely.  The  nonlinear  model  predicted  failure  at 
approximately  1.4  inches.  Actual  failure  took  place  from  2.081  to 
2.114  inches.  This  varied  only  slightly,  compared  to  the  load  at 
failure.  The  load  varied  between  approximately  3600  to  4600  (It  was 
impossible  to  get  an  exact  reading  at  failure  because  samples  were  only 
taken  every  4  seconds).  All  the  curves  were  approximately  the  same 
shape.  The  standard  deviation  for  load  is  356  lbs.  For  displacement 
the  standard  deviation  is  0.012  inches. 

The  failure  occurred  at  so  much  lower  displacement  because 
Capt  Martin’s  data  was  cut  off  at  about  5.6  percent  strain.  The 
computer  assumed  that  at  the  end  of  the  data  failure  had  occurred. 
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Even  with  removing  the  failure  criteria  from  the  program,  it  would  not 
run  past  1.4  inches.  It  was  only  by  extrapolating  the  data  and 
removing  the  failure  criteria  from  the  curves  were  we  able  to  get  the 
results  shown.  To  get  a  better  result  the  Gr/PEEK  should  be  tested 
for  material  properties  up  to  failure  (30%  strain)  and  then  used  in  the 
program. 

The  model  does  closely  follows  the  curve.  The  difference  is 
probably  due  to  the  change  in  material  properties  from  our  specimens 
to  Capt  Martin's  specimens.  In  addition,  the  finite  element  model  is 
naturally  going  to  be  stiffer  than  the  experimental  data.  This  is 
because  one  cannot  have  an  infinite  number  of  elements.  Even  within 
our  own  specimens,  there  is  a  variance  between  the  load  displacement 
curves. 

Plotting  the  same  curve  using  the  notched  specimen  (Figure  5- 
4),  failure  is  predicted  even  earlier.  This  is  caused  by  the  increase  in 
stresses  at  the  point  of  failure.  Otherwise,  the  curve  itself  is  a  very 
close  approximation  to  the  experimental  data. 

This  is  also  the  average  of  seven  specimens  this  time  slightly 
wider  (1.125  ")  with  .1"  radius  notches  on  both  sides  of  the  specimen. 
Displacements  and  loads  are  measured  per  Figure  5- 1  for  the 
experimental  results  and  Figure  5-2  for  the  analytical  results. 
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Displacement  (inches) 

Figure  5-4.  Load  Displacement  Curve  for  Notched  Specimen 


Updating  the  coordinates  is  necessary  at  even  small 
displacements.  This  difference  is  caused  by  the  increased  stress  from 
not  allowing  the  elements  to  displace.  This  is  only  slightly  better 
than  the  linear  program  above  .02"  or  1500  lbs.  This  is  shown  in 
figure  5-5. 
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Figure  5-5.  Load  Displacement  Curves  including  no  updating  of 
Coordinates  for  the  specimen  without  a  discontinuity 

Another  problem  with  the  larger  model  is  the  use  of  triangular 
elements  at  the  point  of  stress  concentration.  This  could  be  part  of 
the  problem  with  the  failure  criteria.  As  the  model  displaces  the  finite 
element  model  also  displaces.  Therefore  a  smaller  model  was  used  to 
show  the  points  of  failure.  This  model  is  discussed  in  Chapter  3  as 
part  of  the  convergence  study. 

B.  Failure  Comparison 

From  Figure  5-6a-d,  the  failure  point  is  at  approximately  1.5 
inches  from  the  tab.  Figure  5-6b  shows  two  failed  specimens  and  two 
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prior  to  testing.  Although  there  is  a  wide  variance  of  ±  1  inch,  this  is 
approximately  where  the  specimen  fails  experimentally.  The 
nonlinear  model  predicted  constant  displacement  along  the  edge. 
This  is  an  obvious  improvement  in  the  linear  model.  The  nonlinear 
model  also  predicted  failure  as  shown. 


Figure  5-6a.  Photograph  of  Failure  of  Plain  Spec  uen 
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Figure  5-6b.  Comparison  of  Specimens 


Failure  occurs  here 


Figure  5-6c.  Finite  Element  Failure  (Small  Model) 


Failure  occurs  at  approximately 


Figure  5-6d.  Experimental  Failures 

From  the  Notched  specimen,  failure  is  predicted  at  slightly  off 
the  center  line  (Figure  5-7a).  This  is  predicted  by  the  nonlinear 
model  (Figure  5- 7b). 


Figure  5- 7a.  Photograph  of  failure  of  Notched  Specimen 


Failure  occurs  here 


Figure  5-7b  Displacement  at  failure  of  Notched  Model 

C.  Displacement.  Strain  Comparisons 

Figure  5-8  shows  the  difference  between  the  strain  (ex )  in  the  x 

direction  and  the  overall  displacement.  Both  sets  of  data  came 
directly  off  of  the  computer  printouts.  Examples  of  the  printouts  are 
shown  in  Appendices  1  and  3. 


Displacement  in  X  Direction  (inches) 


Figure  5-8.  Sc  vs  Displacement:  No  Discontinuity  (Data  as  shown  in 
figure  5-12) 

Figure  5-9  is  a  graph  shows  the  difference  between  the  strain 
(Ey)  in  the  y  direction  and  the  overall  displacement.  This  data  also 

comes  off  the  computer  runs.  One  reason  for  the  differences  in  strain 
for  the  predicted  data  could  be  the  debonding  discussed  earlier  in  this 
section.  This  would  account  for  the  higher  predicted  stain  in  the  x 
direction  and  the  lesser  strain  in  the  y  direction.  Figure  5-10  and 
Figure  5-11  show  much  less  differences.  From  the  experiments  we 
see  much  less  debonding  in  the  notched  specimen. 
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Displacement  in  X  Direction  (in) 

Figure  5-9.  £y  vs  Displacement:  No  Discontinuity  [Data  as  shown  in 
figure  5-12) 

Figures  5-10  and  5-11  show  the  same  data  as  Figures  5-8  and  5- 
9  for  the  notched  specimens. 
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Strain  in  X  Direction  (micro  in/in 


Displacement  in  X  Direction 

Figure  5-10.  Scvs  Displacement:  Notched  Specimen 
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Displacement  in  X  Direction  (in) 

Figure  5-11.  *y  vs  Displacement:  Notched  Specimen 


All  graphs  in  this  section  are  comparing  stresses  and  strains  at  a 

point  1  inch  from  the  tab  (Figure  5-12). 

Comparisons  made  along  this  line 


rosette 

Figure  5-12.  Line  of  Comparison  for  Analytical  Models 

From  Figure  5-13  the  linear  model  is  beginning  to  fall  apart. 
This  is  due  to  the  nonlinearity  of  the  geometry  at  0.5"  displacement. 
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The  nonlinear  program  shows  a  smooth  curve  across  the  specimen. 
The  contour  plotting  also  shows  this. 


Figure  5-13.  Sc  vs  Position  1"  from  tab  for  0.5"  Total  Displacement 
Figure  5-14  also  shows  the  problems  with  the  linear  models. 
Figure  5-15  shows  how  the  shear  strain  goes  to  zero  at  the  line 
of  symmetry  and  at  the  edges  af  the  specimen.  This  is  as  expected 
per  Figure  4-10  we  see  that  shear  strain  is  near  zero  at  the  center 
line.  Figure  5-16  shows  that  shear  strain  does  not  go  to  zero  in  the 
linear  model.  Therefore,  the  nonlinear  program  provides  far  better 
data.  This  appears  to  be  true  for  all  properties. 
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Stress  in  X  Direction 


Figure  5-15.  Yxy  vs  Position  at  0.1",0.5”  and  1"  Displacements 
E.  Contour  Plotting 

Contour  plotting  is  done  using  a  program  written  by  Dr.  Sandhu 
to  show  how  stresses  and  strains  are  varying  with  respect  to  position. 
The  stresses  do  not  vary  per  ply  due  to  the  ±45°  lay-up.  This  is  very 
useful  in  comparing  the  different  types  of  nonlinear  programs.  Figures 
5-17a-d  shows  that  with  the  fiber  orientation  updates  there  is  a 
smoother  transition  across  the  specimen.  Figures  5-17a-b  are  at  a 
displacement  of  1.5".  Figures  5-17c-d  are  at  a  displacement  of  2.5" 
These  also  indicate  a  smoothe.  curve. 


Position  (in) 


Figure  5-16.  Yxy  vs  Position  at  0.5"  Displacement 


Figure  5- 17a.  °x  for  Nonlinear  Program  with  Angle  Updates(1.5") 


lit aao 


Figure  5- 17b.  °x  for  Nonlinear  Program  without  Angle  Updates!  1.5") 


Figure  5- 17c.  cx  for  Nonlinear  Program  with  Angle  Updates(2.5") 
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Figure  5-17d.  ax  for  Nonlinear  Program  without  Angle  Updates(2.5 ") 
Figure  5-18  is  at  the  same  displacement  as  Figure  5-13.  This  is 
also  a  plot  of  strain  in  the  x  direction.  Figure  5-13  goes  across  the 
entire  specimen  where  Figure  5-18  is  only  across  a  quarter  of  the 
specimen. 


Figure  5-18.  ex  for  Nonlinear  Program  with  Angle  Updates(0.5") 

In  addition  to  showing  the  smoother  curve  for  the  plot  of  angle 
updated  data.  Figures  5- 19a  and  b  show  the  strain  concentration  at  the 
edge  upper  left  comer.  This  follows  the  experimental  results,  which 
causes  the  debonding.  These  plots  were  made  at  0.5  inches  which  is 
where  debonding  begins.  Figure  5-20  shows  the  beginning  of 
debonding  at  0.5  inch  displacement  (Look  carefully  at  the  upper 
clamp).  The  strands  of  the  glass  epoxy  can  be  seen  in  the  upper  left 
hand  comer. 


5-19 


-  TOpoo 


Figure  5- 19a.  for  Nonlinear  Program  with  Angle  Updates(0.5") 


Figure  5- 19b.  ty  for  Nonlinear  Program  without  Angle  Updates(0.5 ") 

Figure  5-2 1  shows  the  stress  concentration  slightly  off  the 
center  of  the  circular  discontinuity.  This  rlso  follows  the 
experimental  results  shown  in  Figure  5- 5a. 

Figures  5-22a  and  b  shows  that  the  linear  model  indicates  the 
shear  concentration  at  the  center  of  the  circular  discontinuity.  Since 
it  is  the  shearing  strain  which  causes  the  failure,  the  nonlinear  model 
shows  better  results. 

F.  Displacement  Modeling 

There  is  also  a  plotting  program  to  show  how  the  specimen  is 
displacing  as  the  load  is  applied  by  plotting  the  updated  coordinates. 
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Figure  5-20.  Photograph  at  0.5"  Displacement 


smooth  at  about  one  inch  from  the  tab.  As  the  displacements 
continue,  the  elements  are  less  and  less  rectangular  around  the  one 
inch  point  where  the  failure  occurs. 


Figure  5-23a.  Displacement  for  Nonlinear  Program  with  Angle  Updates 
(0.1") 


Figure  5-23b.  Displacement  for  Nonlinear  Program  with  Angle  Updates 
(0.25") 
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Figure  5-23c.  Displacement  for  Nonlinear  Program  with  Angle  Updates 
(0.5") 


Figure  5-23d.  Displacement  for  Nonlinear  Program  with  Angle 
Updates  (0.75") 


Figure  5-23e.  Displacement  for  Nonlinear  Program  with  Angle  Updates 
(1") 
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Figure  5-23f.  Displacement  for  Nonlinear  Program  with  Angle 
Updates(1.5”) 


Figure  5-23g.  Displacement  for  Nonlinear  Program  with  Angle 
Updates(2") 

A  similar  sequence  of  displacements  is  shown  for  the  notched 
model  in  Figures  24a-c.  Figure  5-24c  is  after  the  specimen  has  failed. 


Figure  5-24a.  Displacements  for  Nonlinear  Program  with  Angle 
Updates(0.5") 


Figure  5-24b.  Displacements  for  Nonlinear  Program  with  Angle 
Updates(l") 


Figure  5-24c.  Displacements  for  Nonlinear  Program  with  Angle 
Update~(2") 
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VI.  Conclusions 


From  the  data  collected,  the  changes  appear  to  provide  the 
program  a  better  method  of  predicting  stresses,  strains,  loads,  and 
failures  for  large  strains.  Updated  Lagrangian  is  essential  for  nonlinear 
materials  at  strains  above  5  percent.  Angle  updating  is  important  at 
strains  above  10  percent.  This  updating  should  in  no  way  hinder  any 
of  the  original  program  capabilities. 

The  use  of  equation  2-9  for  the  calculation  of  the  change  in  fiber 
orientation  is  accurate  even  at  high  strain.  This  is  far  easier  than 
calculating  the  angle  by  examining  the  displaced  elements. 

The  failure  criteria  may  have  several  problems.  For  one,  three 
dimensional  effects  are  not  taken  into  effect.  Another,  the  energy  of 
the  interply  scissoring  is  not  taken  into  account.  A  third  is  that  the 
tab  debonds  at  these  high  strains. 

The  high  elongation  gages  worked  well  even  above  the  their 
stated  maximum.  The  failures  in  the  gages  were  usually  due  to  the 
glue  failing,  the  gages  saturating,  the  gages  peeling  off,  or  the 
connection  breaking.  The  gages  themselves  held  to  failure.  When  the 
specimens  were  tested  the  day  after  applying  the  gages  the  glue  did 
not  fail. 

Recommendations 

By  running  more  experiments  with  high  elongation  gages  this 
data  could  be  even  better.  This  would  provide  better  data  for  the  cubic 
splines  used  by  the  program. 

By  increasing  the  length  of  the  specimen  one  might  avoid  the 
debonding  of  the  tabs.  One  might  also  use  tapered  tabs.  (4)  This 
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Appendix  A:  Program  Inputs,  and  Outputs 

Input  Data  For  Modified  PLSTR2  (small  model): 


PLATE  WITH  NO  DISCONTINUITY 


114 

185 

2  0 

0. 

2 

75 

185  0 

3  1 

1 

11 

0.0 

0.0 

2 

10 

0.0 

0.1 

3 

10 

0.0 

0.2 

4 

10 

0.0 

0.3 

5 

10 

0.0 

0.4 

6 

10 

0.0 

0.5 

7 

1 

0.10000 

0.00000 

8 

0.10000 

0.10000 

12 

0.10000 

0 . 50000 

13 

1 

0.20000 

0 . 00000 

14 

0.20000 

0.10000 

18 

0.20000 

0 . 50000 

19 

1 

.30000 

0  00000 

20 

.30000 

0.10000 

24 

.30000 

0 . 50000 

25 

1 

.41000 

0 . 00000 

26 

.41000 

0.10000 

30 

.41000 

0 . 50000 

31 

1 

.52000 

0 . 00000 

32 

.52000 

0.10000 

36 

.52000 

0 . 50000 

37 

1 

.64000 

0 . 00000 

38 

.64000 

0.10000 

42 

.64000 

0 . 50000 

43 

1 

.77000 

0 . 00000 

44 

.77000 

0.10000 

48 

.77000 

0 . 50000 

49 

1 

.91000 

0  00000 

50 

.91000 

0.10000 

54 

.91000 

0 . 50C00 

55 

1 

1.0600 

0 . 00000 

56 

1  0600 

0.10000 

60 

1  . 06000 

0 . 50000 

61 

1 

1.22000 

0 . 00000 

62 

1.22000 

0.10000 

66 

1.22000 

0 . 50000 

67 

1 

1.39000 

0 . 00000 

68 

1.39000 

0.10000 

72 

1.39000 

0 . 50000 

73 

1 

1.57000 

0 . 00000 

74 

i  57000 

0.10000 

78 

1.57000 

0 . 50000 

79 

1 

1.77000 

0 . 00000 

80 

1.77000 

0.10000 

84 

1.77000 

0 . 50000 

85 

1 

2.00000 

0 . 00000 

86 

2.00000 

0.10000 

90 

2.00000 

0 . 50000 

91 

1 

2.25000 

0 . 00000 

92 

2.25000 

0.10000 

96 

2.25000 

0 . 50000 

97 

1 

2.55000 

0.00000 

98 

2.55000 

0.10000 

0. 
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102 

2.55000 

0.50000 

103 

1 

3.0 

0.0 

104 

3.0 

0.1 

108 

3.0 

0.5 

109 

11 

3.60000 

0.00000 

.025 

no 

10 

3.60000 

0.10000 

.025 

111 

10 

3.6 

.2 

.025 

112 

10 

3.6 

.3 

.025 

113 

10 

3.6 

.4 

.025 

114 

10 

3.60000 

0.50000 

.025 

0.0 

45.  C 

) 

-45.0 

.042 

.042 

.042 

2 

1 

1 

1 

7 

8 

2 

1 

2 

0. 

.042 

15 

11 

12 

6 

5 

1 

-45.0 

.042 

16 

13 

14 

8 

7 

1 

45.0 

.042 

25 

17 

18 

12 

11 

1 

-45.0 

.042 

26 

19 

20 

14 

13 

1 

45.0 

.042 

35 

23 

24 

18 

17 

1 

-45.0 

.042 

36 

25 

26 

20 

19 

1 

45.0 

.042 

45 

29 

30 

24 

23 

1 

-45.0 

.042 

46 

31 

32 

26 

25 

1 

45.0 

.042 

55 

35 

36 

30 

29 

1 

-45.0 

.042 

56 

37 

38 

32 

31 

1 

45.0 

.042 

65 

41 

42 

36 

35 

1 

-45.0 

.042 

66 

43 

44 

38 

37 

1 

45.0 

.042 

75 

47 

48 

42 

41 

1 

-45.0 

.042 

76 

49 

50 

44 

43 

1 

45.0 

.042 

85 

53 

54 

48 

47 

1 

-45.0 

.042 

86 

55 

56 

50 

49 

1 

45.0 

.042 

95 

59 

60 

54 

53 

1 

-45.0 

.042 

96 

61 

62 

55 

55 

1 

45.0 

.042 

105 

65 

66 

60 

59 

1 

-45.0 

.042 

106 

67 

68 

62 

61 

1 

45.0 

.042 

115 

71 

72 

66 

65 

1 

-45.0 

.042 

116 

73 

74 

68 

67 

1 

45.0 

.042 

125 

77 

78 

72 

71 

1 

-45.0 

.042 

126 

79 

80 

74 

73 

1 

45.0 

042 

135 

83 

84 

78 

77 

1 

-45.0 

,042 

136 

85 

86 

80 

79 

1 

45.0 

.042 

145 

89 

90 

84 

83 

1 

-45.0 

.042 

146 

91 

92 

86 

85 

1 

45.0 

042 

155 

95 

96 

90 

89 

1 

-45.0 

042 

156 

97 

98 

92 

91 

1 

45.0 

042 

165 

101 

102 

96 

95 

1 

-45.0 

042 

166 

103 

104 

98 

97 

1 

45.0 

042 

175 

107 

108 

102 

101 

1 

-45.0 

042 

176 

109 

110 

104 

103 

1 

45.0 

042 

185 

113 

114 

108 

107 

1 

-45.0 

042 

17 

16 

15 

20 

15  17 

16 

0. 

0000 

0 

0.0010 

19500. 

0.0020 

0. 

0030 

58500 

0 , 0040 

78100. 

0.0050 

0 

006 

118000 

0.007 

138500. 

0.008 

0. 

009 

180C00 

0.01 

202000. 

0.011 

0 

012 

247500 

0.013 

271000. 

0.014 

0, 

141 

290739 

0.15 

319000. 

0 

0000 

0 

0.0008 

15000. 

0.0016 

0 

0024 

44500 

0.0032 

58500 . 

0  0040 

0. 

0048 

86000 

0.0056 

99000. 

0.0064 
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39000. 

97800. 

159000. 

224500. 

295000. 

30000. 

72500. 

111500. 


would  also  avoid  the  stress  concentration  causing  the  debonding.  This 
would  mean  shaping  the  tabs  in  the  shape  of  the  deformation  shown  in 
Figure  4-9. 
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0  0072 

124000. 

0.0080 

136000. 

0.0088 

148000. 

0 . 0096 

159000. 

0.0104 

170000. 

0.110 

178000. 

0.112 

180500. 

0.0000 

0. 

0.001 

1550. 

0.002 

3100. 

0.003 

4600. 

0.004 

6050. 

0.005 

7450. 

0.006 

8850. 

0.007 

10200. 

0.008 

11500. 

0.009 

12550. 

0.010 

13300. 

0.011 

13900. 

0.012 

14300. 

0.129 

145860. 

0.13 

146000. 

0.0000 

0. 

0.002 

3200. 

0.004 

6400. 

0.006 

9600. 

0.008 

12750. 

0.010 

15500. 

0.012 

18000. 

0.014 

20250 

0 .016 

22200 . 

0.018 

24050. 

0.020 

25730. 

0.022 

27300. 

0.024 

28700. 

0.026 

29900. 

0.028 

30900 . 

0.030 

31800. 

0.032 

32500. 

0.034 

33100. 

0.352 

334000. 

0.36 

336000 . 

0.0000 

0.0 

0.005 

2986.0 

0.0180 

7743.0 

0.034 

9552.0 

0.051 

10420.0 

0.070 

11205.0 

0.108 

12730.0 

0.129 

13517.0 

0.170 

15146.0 

0.209 

16831.0 

0.248 

18530.0 

0.283 

20339.0 

0.318 

21724.0 

0.351 

23186.0 

0.378 

24690.0 

0  0000 

0.305 

0.001 

0.310 

0.002 

0.315 

0.003 

0.3167 

0.004 

0.3165 

0.005 

0.315 

0.006 

0.312 

0.007 

0.3085 

0.008 

0.305 

0.009 

0.3025 

0.010 

0.3 

0.011 

0.298 

0.012 

0  296 

0.013 

0.2943 

0.014 

0.293 

0.500 

0.2929 

0.550 

0.292 

0 . 0000 

0.34 

0 . 0008 

0.35 

0.0016 

0.356 

0.0024 

0.359 

0  0032 

0.3615 

0.0040 

0.3635 

0.0048 

0.365 

0  0056 

0 . 3662 

0.0064 

0.3673 

0  0072 

0.3684 

0 . 0080 

0 . 3695 

0.0088 

0.3705 

0  0096 

0.3715 

0.0104 

0.3725 

0.500 

0.3733 

0  510 

0.3735 

19500000  18750000.  1550000.  1600000.  812500. 

26  26 

26  26 

21  26 

26 

o  oooooo 

0.00 

.001000 

4141.90 

.002000 

8119.04 

003000 

11927.77 

. 004000 

15502.54 

. 005000 

18744.10 

006000 

21988.71 

.007000 

25053.58 

. 008000 

28184.42 

009000 

31303.99 

.010000 

34287.47 

.011000 

37326.15 

012000 

40375.34 

.013000 

43413.02 

.014000 

46419.94 

015000 

49386.09 

.016000 

52321.32 

017000 

55175.65 

018000 

58001.02 

.019000 

60855.01 

. 020000 

63722.05 

021000 

66514.80 

.022000 

69295.50 

. 023000 

72089.79 

0.27573  832174  6  0 

321468  943451  29 

0  OOOOOO 

0.00 

.001000 

4141.90 

. 002000 

8119.04 

.003000 

11927.77 

. 004000 

15502.54 

. 005000 

18744.10 

006000 

21988.71 

.007000 

25053.58 

.008000 

28184.42 

009000 

31303  99 

.010000 

34287.47 

.011000 

37326.15 

012000 

40375.34 

.013000 

43413.02 

.014000 

46419.94 

.015000 

49386.09 

.016000 

52321.32 

.017000 

55175.65 

.018000 

58001.02 

.019000 

60855.01 

.020000 

63722.05 

021000 

66514.80 

022000 

69295.50 

.023000 

72089.79 

0.27573  832174  6  0. 

321468  943451.29 

0  OOOOOO 

0.00 

.001000 

4141.90 

.002000 

8119.04 

003000 

11927  77 

.004000 

15502.54 

.005000 

18744.10 

.006000 

21988  71 

. 007000 

25053.58 

.008000 

28184.42 

009000 

31303  99 

.010000 

34287.47 

.011000 

37326.15 

012000 

40375  34 

.013000 

43413.02 

.014000 

46419.94 

.015000 

49386.09 

.016000 

52321.32 

.017000 

55175.65 

.018000 

58001 .02 

.019000 

60855.01 

.020000 

63722.05 

.021000 

66514.80 

.022000 

69295.50 

.023000 

72089 . 79 
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0  000000 

0.00 

.001000 

4141.90 

.002000 

8119.04 

.003000 

11927.77 

. 004000 

15502.54 

.005000 

18744.10 

006000 

21988.71 

. 007000 

25053.58 

.008000 

28184.42 

.009000 

31303.99 

.010000 

34287.47 

.011000 

37326.15 

.012000 

40375.34 

.013000 

43413.02 

.014000 

46419.94 

.015000 

49386.09 

.016000 

52321.32 

.017000 

55175.65 

.018000 

58001.02 

.019000 

60855.01 

. 020000 

63722 . 0r> 

.021000 

66514.80 

.022000 

69295.50 

.023000 

72089.79 

1.27573  832174.6  0 

.321468  943451.29 

0 . 000000 

0.00 

.005000 

2855.55 

.010000 

4471.86 

.015000 

5442.37 

. 020000 

5973.45 

.025000 

6410.88 

.030000 

6690.98 

. 035000 

6900 . 77 

.040000 

7054  57 

.045000 

7179.76 

. 050000 

7287.89 

.055000 

7396  03 

.060000 

7463.32 

. 065000 

7589.96 

.070000 

7690.74 

.075000 

7749.71 

. 080000 

7895.29 

.085000 

8039.59 

. 090000 

8185.200 

15i743  9430.85  0 

.213487  10676.49 

0  0000000 

1302935 

.001000 

.130293 

.002000 

.129309 

.003000 

.127840 

.004000 

.122488 

. 005000 

.113358 

.006000 

.104879 

.007000 

.097967 

. 008000 

.092782 

.009000 

.087801 

.010000 

.084354 

.011000 

.081168 

.012000 

.078120 

.013000 

.074784 

.014000 

.071385 

.015000 

.067775 

.016000 

.063838 

.017000 

.058719 

.018000 

.054426 

.019000 

.050026 

.020000 

.045578 

.021000 

.041045 

.022000 

.037599 

.023000 

.034035 

27573  0 

.25534  0 

.321468  0 

.17033 

0 . 0000000 

.1302935 

.001000 

.130293 

. 002000 

.129309 

. 003000 

.127840 

. 004000 

.122488 

.005000 

.113358 

. 006000 

104879 

. 007000 

.097967 

. 008000 

.092782 

009000 

.087801 

.010000 

.084354 

.011000 

.081168 

012000 

.078120 

.013000 

.074784 

.014000 

.071385 

015000 

067775 

016000 

.063838 

.017000 

.058719 

01800C 

.054426 

019000 

.050026 

. 020000 

.045578 

021 00C 

.041045 

. 022000 

.037599 

.023000 

.034035 

27573  0. 

25534  0 

321468  0 

17033 

4141900.  4141900.  4141900.  4141900.  571110. 
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Output  Data: 


1PLATE  WITH  NO  DISCONTINUITY 

NUM3ER  OF  NODAL  POINTS -  114 

NUMBER  OF  ELEMENTS -  185 

NUMBER  OF  DIFF.  MATERIALS -  2 

NUMBER  Or  PRESSURE  CARDSS -  0 

X- ACCELERATION -  O.OOOOE+OO 

Y- ACCELERATION -  O.OOOOE+OO 

REFERENCE  TEMPERATURE -  O.OOOOE  +  OO 

MATERIAL  COMBINATION 

**  ISO.-:  ANISO. -2  BOTH-3  **  2 

NUMBER  OF  ISOTROPIC  MATERIALS  0 

TYPE  CF  OUTPUT -  0 

....  ..  .....  NCONTR -  1 

NUMBER  OF  ITERATION -  0 

NONLINEAR  ANALYSIS  NONLIN-1  — 

LINEAR  ANALYSIS  NONLIN-O - 

NCNLIN -  1 

NCPSHN -  2 

MA.X.  NO.  OF  INRETS-MAXINR -  75 

MAXIMUM  ELEM.  OUTPUT  MAXEL -  185 

MA.X .  NC.  CF  INRETS-MAXINR -  75 

MAXIMUM  ELEM.  OUTPUT  MAXEL -  185 

MINIMUM  ELEM.  OUTPUT  MINE- -  0 

DEPTH-WISE  ELEMENTS  MAXDPT -  3 

MDIS^O  FORTE  LOADING - 


I  DISPLACEME 

NT  load:? 

So”” 

cr  lcacing 

MDIS - 

1 

pz:  NT 

TYPE  X 

ORDINATE 

Y  ORDINATE 

X  LOAD  OR  DISPLACEMENT 

Y  LOAD  OR  DISPLACEMENT 

TEMPERATURE 

i 

11 

C.OOCOO 

0.00000 

O.OOOOOOOE+OO 

O.OOOOOCOE+OO 

c .  ccc 

2 

10 

0.00000 

0. 100C0 

O.OOOOOOOE+OO 

O.OOOOOOOE+OC 

o.oo: 

3 

1  c 

0.00000 

0.20000 

O.OOOOOOOE+OO 

O.OOOOOOOE+OC 

c.ccc 

4 

-  A 

C.  00000 

0.30000 

O.OOOOOOOE+OO 

O.OOOOOOOE+OC 

c.oo: 

i 

:  c 

o. ooccc 

0.40000 

O.OOOOCOOE+OO 

O.OOOOOOOE+OC 

C  CC' 

t 

:o 

0.00000 

0.50000 

O.OOOOOCOE+OO 

O.OOCOOOOE-OC 

c.ccc 

- 

i 

0.10000 

0.00000 

O.OOOOOOOE+CO 

O.OCOOOOOE+OC 

0.000 

6 

c 

o. i  coco 

0.10000 

0. OOOCOOOE+CO 

O.OOCCCOOE+OO 

A  AAA 

0  .  V  ^  V 

a 

r\ 

a  a.  a  a  a 

U  .  1  U  w  V 

0.20000 

O.OOOOOOOE+OO 

O.OOOOOCOE+OO 

A  AA+> 
V  •  V  w  -• 

1  3 

c 

C  .10000 

0. 3C00C 

O.OOOOOOOE+OO 

O.OOOOOCCE+CO 

0.000 

1 1 

c 

C.1C00C 

0.40000 

O.OOOCOOCE+OO 

O.OCOCOOOE+OO 

0.000 

1 2 

o 

C.  IC00C 

0 . 5  0  yj  V  O 

0.C00CC00E+CC 

0.0C00000E+C0 

c .  ccc 

3 

1 

0.20000 

c.occcc 

O.COCOOOCE+OO 

O.OCOOOOOE+OC 

c.ccc. 

:  4 

c 

0.2  0000 

0.10000 

C.COCCCOOE+OC 

c.oocooooE-o: 

0  .000 

c 

0.20000 

0.20000 

0.C00CC0CE+00 

O.OOOOOOOE+OO 

c.ccc 

:  e 

c 

C  .20000 

0. 30000 

0.0000CC0E+00 

O.OOOOOOOE+OO 

c.ccc 

: ~ 

r> 

0.20000 

0.40000 

O.OOOOOCOE+OO 

O.OOOOOOOE+OO 

c.ccc 

1  6 

0 

0.2C000 

0.50000 

O.OOOOOOOE+OO 

O.OCOCOOOE+OO 

c.ccc 

1 9> 

l 

0.30000 

0.00000 

O.OOOOOOOE+OO 

O.CCCCCOOE+OO 

c.  c  c . 
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After  all  the  nodes  and  elements  are  shown,  the  program  writes  the  cubic  spline  data: 


BAND  WIDTH 


16 


FOR  THIS  PROGRAM  THE  LOCATION  OF  AA  USED  IS  - 

1 


9411  AND  IN  IA  IS  - 
MATERIAL  1 


STRAIN  0  DEG. (TEN) 
O.OOOOOOOOE+OO 
0 . 1 OOOOOOOE-02 
0.20000000E-02 
0.30000000E-02 
0.4 OOOOOOOE-02 
0 . 50000000E-02 
0 . 60000000E-02 
0 . 7  0000000£-02 
0.80000000E-02 
0 . 90000000E-02 
0 . 1 C0C00CCE-01 
0.  UOOCOOCE-Ol 
0. 120CC0CCE-C1 
0 . 1 3CCO0OCE-O1 
0.14C00C00E-01 
0.141 OOOOOE-Ol 
0 . 15CCC000E-01 


STRESS  0  DEG. (TEN) 
0.00000000E+00 
0 . 19500000E+05 
0 . 39000000E+05 
0.58500000E+05 
0.78100000E+05 
0. 97800000E+05 
0.11800000E+06 
0. 13850000E+06 
0 .  15900000E+06 
0 . 1 8000000E+06 
0.202000CCE+06 
0 . 224  50000E+06 
0.24750000E+06 
0 . 271 OOOOOE+06 
0.29500000E+06 
0 . 29739300E+06 
0. 31900000E+06 


STRAIN  0  DEG. (COM) 
O.OOOOOOOOE+OO 
0 . 8D000000E-03 
0 . 1 6000000E-02 
0 . 24000000E-02 
0 . 32000000E-02 
0.40000000E-02 
0 . 4  8000000E-02 
0 . 56000000E-02 
0 . 64000000E-02 
0.72000000E-02 
0.80000000E-C2 
0 . B8000000E-02 
0 . 96000000E-02 
0. 10400000E-01 
0. 11000000E-01 
0 . 1 120000CE-01 
O.OOOOOOOOE+OO 


1412 


STRESS  0  DEG. (COM) 
O.OOOOOOOuE+OO 
0. 15000000E+05 
0 . 30000000E+05 
0. 44  500000E+05 
0 . 58500000E+05 
0 . 72500000E+05 
0 . 86000000E+05 
0 . 99000000E+05 
0. 11150000E+06 
0 . 12400000E+06 
0 . 1 3600000E+06 
0. 14800000E+06 
0. 15900000E+06 
0 . 1 7000000E+06 
0. 17800000E+06 
0 . 18050000E+06 
O.OOOOOOOOE+OO 


RAIN  9 1 1 E 3 .  ( TEN ) 

STRESS  90DEG. (TEN) 

STRAIN  90DEG . (COM) 

STRESS  90DEG . (COM) 

SHEAR  STRAIN 

SHE 

I . CllCIllCE+CI 

O.OOOOOOOOE+OO 

O.OOOOOOOOE+OO 

O.OOOOOOOOE+CO 

O.OCOCCOCCE-CO 

0.0( 

:.ioo::oooe-c2 

0 . 1 5500000E+04 

0.20000000E-02 

0 . 32000000E+04 

0 . 5O0C00C0E-C2 

0.2! 

o.2:o::oooe-02 

0.31 000000E+04 

0 . 4  OOOCOOOE-02 

0. 64000000E+04 

0. 18000000E-01 

0.71 

o . 3c:::ooce-02 

0 . 4  6000000E+04 

0.60000000E-02 

0 . 96000C00E+D4 

0 . 34  0CC0CCE-C1 

0.9! 

: .4:::::o:e-02 

0. 60500000E+C4 

0.80000000E-02 

0.12750000E+05 

0 . 51 000000E-01 

0 . 1  ( 

: .!::::::cz-02 

0.74500000E+04 

0. lOOOCOOCE-Ol 

0 . 1 5500000E  +  05 

0  - 700000CCE-01 

0.1) 

C -till I II IE-C2 

C . 88S00C00E+C4 

0 . 120COCOOE-C1 

0. 18C00000E+05 

0. 10800000E-00 

o.  i; 

: . -::ci::ce-c2 

0 . 10200CCOE+05 

0. 1400000CE-01 

0 . 20250000E+05 

0 . 12  90000  0E  +  00 

o.i: 

o. &:::::cce-02 

0 . 1 1500000E+C5 

0. 16C00000E-C1 

0 . 22200000E+05 

0. 1700000 0E-00 

0.1! 

: . IE- C2 

0 . 12550000E+05 

0. 180000C0E-C1 

0.24  050000E  +  05 

0.20900000E-0: 

0.1( 

I. 1IIICIICE-C1 

0. 1330000CE+0S 

0 . 20000000E-01 

0 . 257  30000E  +  05 

0.24800000E+00 

0.16 

i . : : :  m:  ce-o: 

0 . 1 39000OCE+05 

0 . 22000000E-01 

0.27  300000E+  05 

0.28300000E-00 

0 . 2  C 

c . :;::::::e-c: 

C . 1 4  300000E  +  05 

0. 24000000E-01 

0 . 287  00000E  +  05 

0 . 31 800000E+00 

0.2) 

I . 12 9III ICE-11 

0 . 1  4  5860C0E  +  05 

0 . 2  6000000E-01 

0 . 29900000E+05 

0.3510000CE+00 

0.2.' 

: . :  oe-oi 

C . 1 4  60000CE  +  05 

0 . 28000000E-C1 

0. 30900000E+05 

0 . 3780C000E+  00 

0.2' 

I.11IICICCE*IC 

O.OOOOOOOOE+CO 

0. 3000C000E-01 

0. 31800000E+05 

0 . OCCCCCOOE^CC 

0.  OC 

O.OOOOOOOOE+OO 

C.32000000E-01 

0 . 32500000E+05 

O.OOOOOOOOE+OO 

0.00 

2.CICII1IIE+CI 

O.OOOOOOOOE+OO 

0 . 34000000E-01 

0. 33100000E+05 

O.OOOOOOOOE+OO 

c.cc 

c .c:ccc:cce*cc 

O.OOOOOOOOE+OO 

0.35200000E-01 

0. 33400000E+05 

O.OOOOOOOOE+OO 

c.o: 

0 . OOOOOOOOE+OO 

O.OOOOOOOOE+OO 

0. 36000000E-01 

0. 33600000E  +  05 

0.00000CCCE+0C 

c.cc 

STRAIN  0  DEG. (TEN) 
O.OOOOOOOOE+OO 
0. 1OCCOOO0E-O2 
0.20000000E-02 
C . 300CC00CE-02 
0. 400C0000E-C2 
0. 5COCCOOOE-02 
0.6000000CE-02 
0 . 7  C00C0C0E-02 
0 . 80000000E-02 
0.90000000E-C2 

.  i  wJJwwJWL 

C.llOOOOOCE-OI 


TEN.  POSSCNS  RATIO 
0.30500000E+00 
0 . 31000000E+00 
0. 31500000E+00 
0 . 31 670000E+00 
0. 31650000E+CC 
0 . 31 500000E+00 
0. 31200000E+00 
0. 30850000E+00 
0 . 30500000E+00 
0.30250000E+00 
0. 30000000E+00 
0 . 2  9800000E  +  00 


STRAIN  90DEG . (COM)  I 
O.OOOOOOOOE+OO 
0 . 80000000E-03 
0. 16000000E-02 
0.24000000E-02 
0. 32000000E-02 
0 . 40000000E-02 
0 . 4  8000000E-02 
0 . 56000000E-02 
0 . 64000000E-02 
0.72000000E-02 
0 . 80000000E-02 
0 . 88000000E-02 


.  POSSONS  RATIO 
0. 34000000E+00 
0. 3S0000C0E+00 
0. 356000C0E+0C 
0. 3S900000E+00 
0. 36150000E+0C 
0. 36350000E+00 
0. 36500000E-00 
C. 36620000E+00 
0. 3673000CE+00 
C . 36840CC0E+C0 
C . 369S0CCCE  +  CD 
0. 3705C0OCE-CC 


A- 6 


0.12000000E-01 
0 . 1 3000000E-01 
0 . 1 4000000E-01 
0 . 14 100000E-01 
0.15000000E-01 


0 . 29600000E+00 
0 . 294  30000E+00 
0 .  29300000E+Q0 
0 . 292  90000E+00 
0.29200000E+00 


0 . 96000000E-02 
0.10400000E-01 
0.11000000E-01 
0.11200000E-01 
O.OOOOOOOOE+OO 


0 .  37 150000E+00 
0. 37250000E+00 
0 . 37  330000E+00 
0 . 37  350000E+00 
O.OOOOOOOOE+OO 


INITIAL  MODULI  OF  ELASTICITY 

E1T=  0 . 19500000E+08  E1C-  0 . 18750000E+08  E2T-  0 . 15500000E+07  E2C-  0 . 16000000E+07  G12» 

MATERIAL  ENERGY-LLT  ENERGY-LLC  ENERGY-TTT  ENERGY-TTC 

1  0 . 20101491E+04  0 . 10398971E+04  0 . 1 1 331 091E+03  0. 75100690E+03 


0.81250000E 

EE 

0 . 538( 


MATERIAL  2 


STRAIN  0  DEG. (TEN) 

STRESS  0  DEG. (TEN) 

STRAIN  0  DEG. (COM) 

STRESS  0  DEG. (COM) 

O.OOOOOOOOE+OO 

O.OOOOOOOOE+OO 

O.OOOOOOOOE+OO 

O.OOOOOOOOE+OO 

0. 10000000E-02 

0. 41419000E+04 

0 . 10000000E-02 

0 . 41 4 19000E+04 

0 . 20000000E-02 

0. 81190400E+04 

0 . 20000000E-02 

0 .811904  00E+04 

0.3C000000E-02 

0. 11 927770E+05 

0 . 30000000E-02 

0.11 927770E+05 

0.40000000E-02 

0 . 15502540E+05 

0 . 40000000E-02 

0.15502540E+05 

0 . 50000000E-02 

0.1B744100E+05 

0 . 50000000E-02 

0 . 18744 100E+05 

0.6000000CE-02 

0.21988710E+05 

0.60000000E-02 

0 . 21 988710E+05 

0.70000000E-02 

0 . 2S053580E+05 

0.70000000E-02 

0.250S3580E+05 

0 . 80000000E-02 

0 . 28184  420E+05 

0.80000000E-02 

0. 2818442 0E+05 

0. 90000000E-02 

0. 31303990E+05 

0 . 90000000E-02 

0. 31 303990E+05 

0. lOOOOOOOE-OI 

0. 34287470E+05 

0. 10000000E-01 

0. 34287470E+0S 

0 . 1 10CC00CE-C1 

0 . 37  326150E+05 

0.11 OOOOOOE-Ol 

0.3732615CE-05 

0 . 1 2000000E-01 

0 . 4  0375340E+05 

0 . 12000000E-01 

0 . 40375  34  0E  +  05 

0 . 1 3000000E-01 

0.4341 3020E+05 

0 . 1 3000000E-01 

0.4341 3020E+05 

0.1400000CE-01 

0 . 4  64 1 994  0E+05 

0. 14003000E-01 

0.464 1 9940E+05 

0 . 1 503000CE-01 

0 . 4  9386090E+05 

0.15000000E-01 

0 . 4  9386090E  +  05 

0 . 1 600000CE-01 

0 . 52321320E+05 

0 . 1 6000000E-01 

0 . 52  321 320E+05 

0.1700000CE-01 

0 . 55175650E+05 

0.17000000E-01 

0 . 55 175650E+05 

0 . 1 800C0C0E-01 

0. 58001020E+05 

0 . 18000000E-01 

0 . 58001020E+C5 

0 . 1 900000CE-01 

0 . 60855010E+05 

0 . 1 9000000E-01 

0. 60855010E+05 

0 . 2  CCCC00CE-01 

0.63722050E+05 

0 . 20000000E-01 

0 . 63722050E+05 

0 . 2 lOOCCOOE-Ol 

0 . 66514800E+05 

0.21000000E-01 

0 . 66514  800E+05 

0.220000COE-01 

0 . 69295500E- 05 

0.22000000E-01 

0 . 692  9S500E+05 

0 . 2  300CC00E-01 

0.72089790E-C5 

0 . 2  3000000E-01 

0 . 72  0897  90E+  OS 

0.27S73CCCE-00 

0 . 8321 74  60E»  06 

0 . 2757  3000E+  00 

0 . 832 1 74  6CE+06 

C.321468CCE-00 

0.94  34  5129E  +  06 

0. 32146800E+00 

0 . 94  34512  9E  +  06 

' 

:  9: -EC* .  (TEN) 

STRESS  90DEG. (TEN) 

STRAIN  90DEG . (COM) 

STRESS  90DEG .  (COM) 

SHEAR  STRAIN 

SHE 

£ 

.  ZZZZZZZZZ-ZZ 

O.OOOOOOOOE+OO 

O.OOOOOOOOE+OO 

O.OOOOOOOOE+OO 

O.OOOCOOOOE-CO 

0  OC 

z 

-C2 

0.4141 9000E  +  04 

0. 10000000E-02 

0. 41419000E+04 

0 . 5000CCC  0E-C2 

0.26 

V 

.2:::::::e-02 

0 . 8 1  904  00E+  04 

0.20000000E-02 

0. 81190400E+04 

0.100000CCE-01 

0.4' 

z 

. 3C: ZZZZZZ-Z2 

0.11 927770E*05 

0. 30000000E-C2 

0.11 927770E*  05 

0.1500CCOOE-01 

0.5' 

z 

0. 15502540E+05 

0 . 4  C00000CE-02 

0 . 1 55C254  0E  +  C5 

0 . 20000000E-C1 

0.5? 

z 

.tz::z:zzz-z2 

0.1874  4 100E+05 

0 . SOCCOOOOE-02 

0. 18744100E+C5 

0.25COOOCOE-01 

0.6' 

z 

.  6:::::::e-02 

C . 2 1 9887 iCE+05 

0.60CC00C0E-02 

0. 2198871 0E+C5 

0 . 30C00C00E-01 

0.  6( 

z 

.  "ZZZZZZZZ -Z2 

C.25053580E*05 

0.70000000E-02 

0 . 25053580E+  05 

0. 35000000E-:: 

C.6! 

z 

.SZZZZZZZZ-Z 2 

0.281 84420E+05 

0 . 80000000E-02 

0.28184420E+05 

0 . 40000CCCE-C 1 

0. 7( 

c 

.  9ZZZZZZZZ-Z2 

0.31 303990E+05 

0 . 90000000E-02 

0. 31303990E+05 

0 . 450000CCE-C 1 

0.7) 

c 

. :  zzzzzzzz-zi 

0. 34287470E-0S 

0 . 10000000E-01 

0. 34287470E+0S 

0 . 50000000E-C 1 

c . 

z 

.  IICCCCCCE-Cl 

0 . 37326150E+05 

0. 11000:70E-01 

0. 37326150E+05 

0 . 550000CCE-01 

c.  7; 

z 

.  i2c::::ce-ci 

0. 4037534 0E+05 

0 . 12000000E-01 

0.40375340E+05 

0. 600000CCE-01 

0.7' 

z 

.  13CCC000E-01 

0.4341 3020E+05 

0 . 1 3000000E-C1 

0.434 1 3020E+  05 

0. 650000C0E-C1 

c.7; 

c 

. :  zzzzzzzz-zi 

0. 46419940E+05 

0 . 14000000E-01 

0. 46419940E+05 

0.7000000CE-01 

C.7( 

c 

.  ISOOOOOOE-Ol 

0 . 49386090E+05 

0 . 1 5000000E-01 

0.49386090E+05 

0.75000000E-C1 

0. 71 

o 

.  :6::oocce-ci 

0.52321320E+C5 

0 . 1 6000000E-01 

0 . 52  32 1 320E  +  05 

o.eocococoE-:: 

0 . 7  ( 

c 

.  17CCC000E-01 

0 . 551 75650E+05 

0.17000000E-0! 

0.55175650E+05 

0 . 85000000E-01 

0. 8( 

c 

.  :s::::ocE-oi 

0 . 58001020E+C5 

0 . 1 8000000E-01 

0 . 58001020E+  05 

0 . 90000000E-C1 

0.8) 

z 

. :  9c::::ce-oi 

0 . 6085501 0E+05 

0 . 1 9000000E-01 

0. 60855010E+05 

0. 15174300E‘CC 

0.9' 

z 

.2o:::c:oe-ci 

0.63722050E+05 

0.20000000E-01 

0.63722050E-05 

0.21348700E*:: 

0 . 1  ( 

c 

. ; :  :c:ccoe-ci 

C.66514800E-05 

0 . 21 OOOOOOE-Ol 

0. 66514800E*05 

O.OOCOCCCCE-OC 

0. 0( 

c 

.22o:::cce-ci 

0 . 6  9295500E  +  05 

0 . 22000000E-01 

0. 69295S00E-C5 

o.cccccccce-c: 

0. 0( 

c 

.  2  3 ZZZZZZZ-ZZ 

0 . 7  20897  90E+  05 

0.23000000E-01 

0.72089790E-05 

o.cococ:::e-:: 

0.0' 

c 

.  2 " 5 7 3CC2E  +  00 

0.83217  4  60E  +  06 

0.27573000E*C0 

0. 83217460E+C6 

o.occcccc:* 

c.  o; 

0 . 32 1 46800E+00  0. 94345129E+06  0. 321468C0E+00  0 . 94 345129E+ 06  0 . COOOOOOOE+CO  O.Ot 


STRAIN  0  DEG. (TEN) 
O.OOOOOOOOE+OO 
0 . 10000000E-02 
0 . 20000000E-02 
0 . 30000000E-02 
0 . 4  0000000E-02 
0.50000000E-02 
0.60000000E-02 
0 . 7  0000000E-02 
0.80000000E-02 
0 . 90000000E-02 
0 . 10000000E-01 
o.nooooooE-oi 
0.12000UUOL-01 
0 . 1 3000000E-01 
0 . 1 4000000E-01 
0. 15000000E-01 
0 . 1 6000000E-01 
o .  nooooooE-oi 
0 . 1 8C0000CE-01 
0.19000000E-01 
0.20000000E-01 
0 . 2IC0000CE-01 
C.2200CC0CE-01 
0 . 2  3000CGOE-01 
C.27573CC0E+00 
0.180CCC03E-C1 
0.:9CC00C0E-C1 


TEN.  POSSONS  RATIO 
0. 13029350E+00 
0.13029300E+00 
0. > 2930900E+00 
0.12784000E+00 
0.12248800E+00 
0 .  11335800E+00 
0. 10487900E+00 
0.97967000E-01 
0.92782000E-01 
0. 87801000E-01 
0.84354000E-01 
0 . 811 68000E-01 
0 . 78120000E-01 
0.74784000E-01 
0.71385000E-01 
0.67775000E-01 
0 . 63838000E-01 
0.58719000E-01 
0.54428000E-01 
0.50026000E-01 
0.45578000E-01 
0.4104S000E-01 
0.37S99COOE-01 
0 . 34035000E-01 
0 . 2SS34  OOOE+OO 
0.54426000E-01 
0 . 5C026000E-01 


STRAIN  90DEG. (COM) 
0.00000000E+00 
0. 10000000E-02 
0.20000000E-02 
0 . 30000000E-02 
0. 40000000E-02 
0 . 50000000E-02 
0.60000000E-02 
0.70000000E-02 
0 . S0000000E-02 
0 . 90000000E-02 
0. 10000000E-11 
0.11000000E-01 
0.12000000E-01 
0 . 1 3000000E-01 
0.14000000E-01 
0 . 15000000E-01 
0 . 1 6000000E-01 
0 . 17000000E-01 
0 . 18000000E-01 
0 . 1 9000000E-01 
0 . 20000000E-01 
0.21000000E-01 
0.22000000E-01 
0.23000000E-01 
0.27573000E+00 
0 . 1 8000000E-01 
0 . 1 9000000E-01 


COM.  POSSONS  RATIO 

0 . 13029350E+00 
0 . 1 3029300E+00 
0 . 12  930900E+00 
0.12784000E+00 
0.12248800E+00 
0. 11335800E+00 
0. 10487900E+00 
0. 97967000E-01 
0.92782000E-01 
0. 8780100 0E-01 
0 . 84  354000E-01 
0 . 81 1 68000E-01 
0 . 78120000E-01 
0 . 74  784000E-01 
0.71 385000E-01 
0.67775000E-01 
0 . 63838000E-01 
0 . 58719000E-01 
0. 54  426000E-01 
0.5002 6000E- 01 
0 . 45578000E-01 
0.41 045000E-01 
0 . 37599000E-01 
0 . 34  035000E-01 
0 . 25534000E+00 
0 . 54  426000E-01 
0 . 50026000E-01 


A-8 


0 . 22000000E-01 
0.23000000E-01 
0.27S73000E+00 
0 . 3214  6800E-I-00 


0.37599000E-01 
0. 34035000E-01 
0.25534000E+00 
0.17033000E+00 


0.22000000E-01 
0.23000000E-01 
0.27573000E+00 
0 . 3214  6800E+00 


0.37599000E-01 

0.34035000E-01 

0.25534000E-00 

0.17033000E*00 


The  program  now  reads  the  Initial  values  of  the  material  properties: 


INITIAL  MODULI  OF  ELASTICITY 

E1T=  0 . 41419000E+07  E1C-  0 . 41419000E+07  E2T-  0 . 4 14 19000E+ 07  E2C-  0 . 4 14 19000E+07  G12  =  0.571110001 

MATERIAL  ENERGY-LLT  ENERGY-LLC  ENERGY-TTT  ENERGY-TTC  El 

2  0.11 639355E+06  0.11639355E+06  0. 11639355E+06  0 . 11 639355E+06  0.11395671 


he  Iterations  are  now  listed  and  the  Increment  given 


ITERATION 

CONTROLS 

ITERATION 

2  CHECK1 

O.OOOOOOOOE+OO 

CHECK2 

ITERATION- 

CONTROLS 

ITERATION 

3  CHECK1 

0 . 12406354E+02 

CHECK2 

ITERATION 

CONTROLS 

ITERATION 

4  CHECK1 

0.13138208E+02 

CHECK2 

ITERATION 

CONTROLS 

ITERATION 

5  CKECK1 

0 . 1 3201 918E-*  02 

CHECK2 

in  the  far  right  hand  column: 


0 . 124  06354E+02 

DECHK 

0.10000000E+01 

SCALE 

o 

o 

o 

o 

o 

0.131 38208E+02 

DECHK 

0 . 55704235E-01 

SCALE 

0.1000C 

0.1320191 BE+02 

DECHK 

0.48258377E-02 

SCALE 

0.1000C 

0 . 1 3213167E+02 

DECHK 

0 . 851 31 919E-03 

SCALE 

0.1000C 

The  program  now  writes  the  new  angles  after  the  first  iteration: 


THE  ANGLES  HAVE  BEEN  CHANGED 

1  C.OOOOOOOOE+OO  2  0 . 44  99"7  4  87E+02 


5 

C . 44991883E+02 

6 

-0 .  44  991883E+02 

9 

-C.4498C874E+02 

10 

O.OOOOOOOOE+OO 

13 

C.OOOOOOOOE+OO 

14 

0 . 44  929219E+02 

17 

-0.4  4  992181E+C2 

18 

0.4497311 1E+02 

21 

-C.4493S856E+02 

22 

0 . 4  4888812E+02 

25 

-0.44  900062E+02 

26 

0.44985753E+02 

29 

-0.4  4  952023E+02 

30 

0 . 4  4  914241E+02 

33 

-C . 4  4  91 6338E+02 

34 

0 .  44  9894S4E+02 

37 

-C.44979148E+02 

38 

0.44  94  3508E+02 

41 

-0.4493987CE+02 

42 

0 . 44  980161E+02 

45 

-C.44997305E+02 

46 

0.44984222E+02 

43 

-0.44970504E+C2 

50 

0 . 44  986742E  +  02 

53 

-C . 44996012E+02 

54 

0 . 44  999032E+02 

57 

-0 . 44997119E-02 

58 

0. 44990347E+02 

61 

-C .  4  4  9  92033E*  02 

62 

0 . 44  996933E+02 

€5 

-0.44  9999*4E-02 

66 

0 . 44  9884  4  7E  +  02 

69 

-0 . 449'7245E-C2 

70 

0 . 44  980688E- 02 

7  2 

-  0 . 44  36594  6E-C2 

74 

0. 44998159E-02 

7  7 

-0.44954168E+02 

78 

0. 44983482E  +  02 

SI 

-:.449*9781E-C  2 

82 

0.44  985037E  +  02 

6: 

-0 . 44595490E-12 

86 

0 . 4  4  997  658E  +  02 

6  9 

-0 . 44392489E-C 2 

90 

0.44987421E+02 

93 

-  C  .  4  <  9  B  6  9  32  E>  C2 

94 

0 . 4  4  993927E+  02 

0  - 

-..44 7  7  9  C  5 . Z  +  22 

98 

0.44997130E+02 

C  1 

.4493313 9E *02 

102 

0.44933728E-02 

Z  5 

-0 .449955183-02 

106 

0.44999980E+02 

:  9 

- t . 4499993. 5  E-C2 

110 

0.44999789E-02 

-0 . 4 4999316E-02 

114 

0.44998645E-02 

i  ~ 

-0 . 44 ;?9375£.02 

118 

0.44998312E-02 

;  i 

-0  .  44  a^SOlE-OO 

122 

0 . 44998163E-02 

2  5 

-0 . 44939531E-02 

126 

0.44999287E-02 

:  ? 

-0 . 44 998066E-02 

1 30 

0.44997505E-02 

-t . 4499'784E*C2 

134 

C.44999012E-C2 

3  ~ 

-2 . 4499965  2Z+Z2 

138 

0.44  999044E-02 

tm  ^ 

-  3 . 4  4  998681 1- C2 

142 

0.44998677E-02 

1  Z 

*.448931025*02 

146 

0.44  999  988E  +  C2 

A  9 

-2 . 44999936E-C2 

150 

0 . 4  4  999825E- 02 

5  2 

-1 .449396*65-02 

154 

0.44999582E-C2 

C  - 

-1 .44333=971-02 

158 

0.44999730E-C2 

6: 

-1 .449396*45-02 

i  €  2 

0.44999769E-02 

Si 

-0 .449999*85-02 

166 

0.44  999954E-C2 

5  ~ 

-2.44  999S97E>32 

158 

0.44999730E-02 

£  \ 

-1.44 9996*4E-02 

162 

0.4499976  9E-02 

1 5 

-0 .449999785-02 

166 

0.44999954E-02 

6  9 

-0.449996685-02 

1 7  C 

0.44  999805E-02 

?  3 

-0 . 44999*81 E- 02 

174 

0.44999814E-02 

7  7 

-0 . 450000COE-C2 

178 

0. 44999998E-02 

e; 

83 

-0.449999905-02 

-0 . 44339330E-C2 

182 

0. 44999973E-02 

3  -0.44997487E+02  4  0 . 00000000E+00 


7 

O.OOOOOOOOE+OO 

8 

0 . 44  980874E  +  02 

11 

O. 44942155E+02 

12 

-0 . 44  942 155E+  02 

15 

-0.44  92921 9E+02 

16 

0 . 44  9921 81E+02 

19 

-0 . 44  97  3111E+02 

20 

0 . 44  938856E+02 

23 

-0 . 44  888812E+02 

24 

0 . 44  900062E+02 

27 

-0 . 44  985753E+02 

28 

0 . 44  952023E+0? 

31 

-0 . 44  91424 1E+02 

32 

0 . 44  91 6338E+02 

35 

-0. 44  989454E+02 

36 

0 . 44  97  91 48E+02 

39 

-0. 4  4  943508E+02 

40 

0 . 44  939870E+02 

43 

-0 . 44  9801 61E+02 

44 

0.44997305E+02 

47 

-0.44984222E+02 

48 

0.44970504E+02 

51 

-0.44986742E+02 

52 

0 . 4499601 2E+02 

55 

-0.44999032E+02 

56 

0 . 44  99711 9E+02 

59 

-0.44990347E+02 

60 

0.44992033E+02 

63 

-0 . 44  996933E  +  02 

64 

0 . 44  99  9974E+02 

67 

-O . 44  9884  47E  +  02 

68 

0 . 44  97  724  5E  +  02 

71 

-0.44980688E+02 

72 

0.4498994  6E  +  02 

75 

-0 . 44  998159E+02 

76 

0.449941 68E+02 

79 

-0.44983482E+02 

80 

0.44979781E+02 

83 

-0 . 44  985037E+02 

84 

0.44  9954  90E  +  02 

87 

-O. 4499765 8E +02 

88 

0 . 4  4  9924  89E+02 

91 

-0.44987421E+02 

92 

0 . 4  4  98  6932E+02 

95 

-0.44993927E+02 

96 

0 . 4  4  999051E  +  02 

99 

-0.4  4  997130E+02 

100 

0 . 44  9951 39E+02 

103 

-0.44993728E-C2 

104 

0 . 44  99551 8E+C2 

107 

-0.4  4  99998CE+  02 

108 

0 . 4  4  99  9935E  +  02 

111 

-0.44999789E+02 

112 

0.44999316E+02 

115 

-0.44998645E+02 

116 

0.44999375E+02 

119 

-0 . 4  4  998312E  +  02 

120 

0.44997801E+C2 

123 

-0.44998163E+02 

124 

0.44999531E-02 

127 

-0.44999287E+02 

128 

0.44998086E-02 

131 

-O. 4 4 997 505E+ 02 

132 

0.44997784E-C2 

135 

-0.44999012E+02 

136 

0.44999652E+C2 

139 

-C.44999044E+02 

140 

0.44998681E+02 

143 

-0.44998677E+02 

144 

0. 449991 C2E+02 

147 

-0 . 44  9999B8E*  02 

148 

0.44  99  9936E+  02 

151 

-O. 44999825E+02 

15.' 

0. 4499967 6E+CL 

155 

-0.44999582E+02 

156 

0.44999897E+02 

159 

-0 . 4  4  9997  30E  +  C2 

160 

0.44999674E+02 

163 

-0 . 4  4  99976  9E  +  C2 

164 

0 . 44  99  9978E- 02 

167 

-C . 44999954E+02 

168 

0.44999868E+C2 

159 

-0.4499973CE+02 

1  60 

0.44999674E+C? 

163 

-0.44999769E+02 

164 

0.44999978E+02 

167 

-O. 44999954E+02 

168 

0.44999868E+02 

171 

-O. 44999805E+02 

172 

0.44999781E-02 

175 

-O. 44999614E-02 

176 

0. 4500CCC0E+02 

179 

-O. 44999998E-02 

180 

0. 44999990E+02 

183 

-0.44999973E+02 

184 

0.44  999950E  +  02 
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e  program  now  writes  the  dispalcements  for  the  Increment: 


LOAD  INCREMENT  1 
DISPLACEMENTS 

N.P.  UX  Of  N .  F .  UX  UY  N.P.  UX 

i  0 . OCOCOOOOE+OO  O.OOOOOOOOE  +  OO  2  O.OOOOOOOOE+OO  -0. 35198089E-03  3  0 . OOOOOOOCE+OO  -0.( 

4  O.OOOOOOOOE+OO  -0 . 10315992E-02  5  O.OOOOOOOOE+OO  -0 . 13142859E-02  6  0 . OOOOOOCOE+OC  -0.1 

7  0. 36206947E-03  O.OOOOOOOOE+OO  8  0 . 361 84309E-03  -0. 35467879E-03  9  0. 35996848E-03  -0.7 

10  0. 35178756E-03  -0. 10408004E-02  11  0 . 35850726E-03  -0. 13791095E-02  12  0 . 53584 429E-03  -0.1 

13  0 . 83956609E-03  0 . 00C00000E+00  14  0 . 84004871E-03  -0. 36403318E-03  15  0 . 842 15591E-03  -0.7 

16  0.86446990E-03  -0 . 110391 31E-02  17  0 . 974 80282E-03  -0. 15624221E-02  18  0 . 121 9094 0E-02  -0.2 

19  0. 13396937E-02  O.OOOOOOOOE+OO  20  0 . 1 34 57754E-02  -0 . 3871 7426E-03  21  0 . 137 74 393E-02  -0.7 

22  0. 14770619E-02  -0 . 12591 805E-02  23  0 . 1 6656396E-02  -0. 18034098E-02  24  0 . 185612 92E-C2  -0.2 

25  0 . 1 94  4  5875E-02  O.OOOOOOOOE+OO  26  0 . 1 96 94 84 4E-02  -0 . 4 47 93927E-03  27  0.20572164E-02  -0.5 

28  C. 22110034 E-02  -0 . 1472017 9E-02  29  0 . 237 69577E-02  -0.20059500E-02  30  0.25378994E-02  -0.2 

31  0.26465411E-02  O.OOOOOOOOE+OO  32  0 . 2 6881159E-02  -0 . 54125260E-03  33  0 . 27  974 61 5E-02  -0.1 

34  C.29317C96E-02  -0 . 164 4 6009E-02  35  0. 30712892E-02  -0 . 21745587E-02  36  0 . 321 1 8812E-C2  -0.2 

37  C . 35262694E-02  O.OOOOOOOOE+OO  38  0 . 35561231E-02  -0 . 62294 459E-03  39  0 . 362 8332 3E-02  -0.1 

40  0.37231587E-02  -0 . 1 7856706E-02  41  0 . 38306582E-02  -0 . 23257020E-02  42  0 . 394 6337CE-02  -0.2 

43  0.45358382E-02  O.OOOOOOOOE+OO  44  0 . 454 18125E-02  -0 . 65204554E-0">  45  0 . 45656780E-02  -0.1 

46  0 . 46109380E-02  -0 . 1 877 1 94 6E-02  47  0 . 4 67493B2E-02  -0 . 24 4 1 9199E-02  48  0 . 4 754 12C2E-02  -0.2 

49  C.56150797E-02  O.OOOOOOOOE+OO  50  0 . 561 22027E-02  -0 . 64 6882 38E-03  51  0 . 5607 391 8E-02  -0.1 

52  0.56099956E-02  -0 . 1906335 3E-02  53  0 . 562 65747E-02  -0 . 250121 64E-02  54  0. 56594425E-02  -0.1 

55  0 . 67458005E-02  O.OC'OOOOOOE+OC  56  0 . 67  4 1061  IE-02  -0 . 631 37  052E-03  57  0 . 672 784 92E-C2  -0.1 

58  0.6^0992341-02  -0 . 1 8860288E-02  59  0 . 66941553E-02  -0 . 24 999308E-02  60  0 . 6687 3201 E-02  -0.2 

61  C  .  ->92  41  721E-C2  0 . 000000C0E  + 00  62  0 . 791  981 4 9E-02  -0. 61504982E-03  63  0 . 790 704  33E-02  -0.1 

64  0. 788'0539E-02  -C . 1 84 691 08E-02  65  0 . 78627326E-02  -0 . 24 632287E-C2  66  0 . 7838990 3E-02  -0.2 

67  C. 51  32  6  94  9 E-02  0 . OOOOOOOOE  +  CC  68  0 . 914 97429E-02  -0 . 60258552E-03  69  0 . 91 4 1094 6E-02  -0.1 

7 -  0 . 912734  30E-C2  -C . 1 81 32 1 06E-02  71  0 . 91 0  94 397E-02  -0 . 24237 952E-02  72  0 . 9088' 97  4E-02  -0.2 

'3  0. 10439456E-C1  C . OOOOCCOOE+OO  74  0. 104 38221E-01  -0 . 5963574 0E-03  75  0 . 104 34 58 3E-01  -0.1 

76  C.10428676E-CI  -0 . 17947807E-C2  77  0 . 104 2C528E-01  -0 . 2 3991 8C6E-02  78  0 . 104 09931E-C1  -0.2 

0 . 1 1 865S28E-C1  O.OOOOOOOOE+OO  80  0 . 1 1 86587 9E-01  -0 . 59564859E-03  81  0 . 1 1 8 65935E-01  -0.1 

S2  C.118657;:e-0I  -0. 179C3256E-02  83  0. 1 1864762E-01  -0.23905981E-02  84  0. 1 1 86254 4E-01  -0.2 

8:  0. 13510346 E-CT  0 . 00C000C0E  +  C0  86  0 . 1 35 10824E-01  -0 . 59779583E-03  87  0 . 1 35 1 2 1 4 5E-01  -0.1 

8-  0. 1  3:  1  4004  E- 01  -0 . 1 7 94 039 6E-02  39  0 . 1 351 5987E-01  -0 . 2392 7986E-02  90  0 . 1 35 1 7 6 3 9E-0i  -0.2 

9  1  0 . 1  5  30  3372  E-01  O.OOOOOOOOE  +  OO  92  0 . 1530367BE-01  -0 . 5997 9029E-03  93  0 . 1 5304 557E-01  -0.1 

94  0 . 15303905E-01  -0 . 1 7 98877 OE-02  95  0  152.07578E-01  -0 . 23580410E-02  96  0. 153094  34E-01  -0.2 

5'  0 . 1'45S127E-01  C  .  OCOOOOOOE  +  OO  98  0 . 1 74 581 94E-01  -C . 60054 67 3E-0 3  99  0 . 17 4 58 3 97E-0I  -0.1 

0.1'  4  56'42E-01  -0 . 1 8012669E-02  101  0 . 1 74  5922  6E-01  -0 . 2401 31 66E-02  102  h  .  1  74  598  34E-01  -0.2 

0  .2  :  6  903  855-01  0.  COOOOCOOE-OO  104  0.20690564E-01  -0 . 6005091 0E-03  105  0 . 20  6  90507E-01  -0.1 

2  .::  6  9:  4  3'E-01  -C.18014'0  3E-C2  1C7  C  .  206  9038GE-C1  -0 . 24  01894  6E-C2  108  0 . 206  90368E-01  -0.2 

229  : .IIOOOOOOE-Ol  O.OOOOOOOOE+OO  HO  0.25CO:OOCE-01  -C.60C3835SE-C3  111  0 . 25000C00E-C 1  -c.i 
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The  program  now  prints  the  stresses,  strains,  and  failure  criteria  numbers: 


1 


STRESSES  /  STRAINS  /  ENERGY  CONTRIBUTIONS 


£W. 

X 

Y 

-XX 

-YY 

-XY 

-LL  -TT  -LT  ENER 

.  LEVi 

1 

0.050 

0.050 

0 .  1272E+05 

0. 3620E-02 

-0. 1235E+05 

-0.3533E-02 

-0. 9210E+01 
-0.1462E-04 

0. 1272E+05  -0. 1235E+05  -0.9210E«01 

0 . 3620E-02  -0 . 3533E-02  -0.1462E-04 

0 . 2265E-03  0.2161E-03  0.5908E-07 

0.441 

2 

C.050 

0.050 

0 . 4714E+04 

0 . 3620E-02 

-0.3627E+04 
-0. 3533E-02 

-0 . 4S40E+03 
-0.1462E-04 

0.10U8E+04  0 . 7982E+02  0.4170E+04 
0.a044E-04  0.3582E-04  0.7153E-02 

0. 1234E-04  0.8760E-05  0.2835E-02 

0.28! 

3 

0.050 

0.C50 

0 . 4582E+04 

0. 3620E-02 

-0. 3759E+04 
-0. 3533E-02 

0. 3160E+03 

-0.14  62E-04 

0.7276E+03  0.9565E+02  -0.4170E+04 

0. 3582E-04  0.5044E-04  -0.7153E-02 

0. 6223E-05  0.1737E-04  0.2835E-02 

0.28! 

4 

0.050 

0.149 

0 . 1270E+05 

0 . 3609E-02 

-0 . 1217E+05 

-C.3484E-02 

-0.2974E+02 

-0.4722E-04 

0.1270E+05  -0 . 1217E+05  -0.2974E+02 

0. 3609E-02  -0 . 3484E-02  -0.4722E-04 

0.2252E-03  0.2102E-03  0.6163E-06 

0.43! 

5 

0.050 

0.149 

0 . 504  9E+04 

0 . 3609E-02 

-0.3231E+04 

-0.3484E-02 

-0 . 8067E+03 
-0.4722E-04 

0.1716E+04  0 . 1021E+03  0.4140E+04 
0.8639E-04  0.3917E-04  0.7093E-02 

0. 3620E-04  0.1047E-04  0.2789E-02 

0.28! 

6 

C.050 

0.149 

0 . 4  622Et04 

0 . 3609E-02 

-0 . 3658E+04 
-0.3484E-02 

0. 3286E+03 
-0 . 4722E-04 

0 . 8104E+03  0 . 1533E+03  -0.4140E+C4 

0 . 3917E-04  0 . 8639E-04  -0.7093E-02 

0.7440E-05  0 . 5095E-04  0.2789E-02 

0.28' 

7 

C.05C 

0.249 

0.12S7E+05 

0 . 3559E-02 

-0.1168E+0" 

-0. 3345E-02 

-0 . 7009E+02 

-0 . 1113E-03 

0 . 1257E+05  -0 . 11 68E+05  -0.7009E+C2 

0 . 3559E-02  -0.3345E-02  -0.1113E-03 

0.2191E-03  0 . 1942E-03  0.3422E-05 

0.41! 

e 

”'”5" 

0.249 

0.5730E+04 

0 . 3559E-02 

-0 . 2358E+04 

-0.3345E-02 

-0. 1529£*04 
-0 . 1113E-03 

0. 3215E+04  0.1571E+03  0.4044E-*04 

0 . 1 624E-03  0.5116E-04  0.6904E-02 

0.1280E-03  0.1787E-04  0.2649E-02 

0.27! 

9 

=.=== 

0.249 

C.4724E*04 

0 . 3559E-02 

-0 . 3365E+04 

-0. 3345E-02 

0. 4023E+03 
-0 . 1113E-03 

0 . 1082E+04  0.2775E+03  -0.4044E-.04 
0.5116E-04  0 . 1624E-03  -0.6904E-02 

0. 1269E-04  0 . 1801E-03  0.2649E-02 

0.28' 

h,Ch 

C.349 

0.1266E-0S 

0. 3551E-02 

-0 . 1077E+05 
-0. 3105E-02 

-0 . 2119E+03 
-0 . 3365E-03 

0 . 1266E+05  -0 . 1077E+05  -0.2119E+03 

0. 3551E-02  -0 . 31 05E-02  -0.3365E-C3 

0 . 2182E-03  0 . 1 67BE-03  0.3129E-04 

0.41'. 

“  “ 

C  .  C5C 

0.349 

0 . 7912E+04 

0 . 3551E-02 

0 . 7754E+02 

-0 . 3105E-02 

-0 . 3722E+04 

-0 . 3365E-03 

0.7717E+04  0. 2724E+03  0.3927E+C4 

0 . 3915E-03  0 . 54  98E-04  0.6656E-C2 

0.7434E-03  0.2064E-04  0.2469E-02 

0.32: 

-  £. 

~  *  w  *"  w 

0.349 

0 . 4  869E+04 

0 . 3551E-02 

-0 . 2966E+04 

-0.3105E-02 

0.3149E+03 
-0 . 3365E-03 

0 . 1266E+04  0. 6365E+03  -0.3917E-04 

0.5498E-04  0.3915E-03  -0.6656E-02 

0.1466E-04  0 . 1 04  7E-02  0.2469E-C2 

0.35! 

CD 

4T 

o 

0 . 1593E-05 

0.4472E-02 

-0 . 1 1 10E+05 

-0 . 3307E-02 

-0 .2593E+03 

-0.4118E-03 

0.1593E+05  -0.11 10E+05  -0.2593E-03 

0.4472E-02  -0.3307E-02  -0.4118E-03 

0 . 34 10E-03  0 . 1899E-03  0.4685E-C4 

0.57'. 

0.448 

0. 1279E+05 
0.4472E-02 

0 . 3836E+04 

-0.3307E-02 

-0 . 7350E+04 
-0.4118E-03 

0 . 1566E+05  0 . 9642E+03  0.4478E+04 

0.7B81E-03  0 . 37  64E-03  0.7779E-02 

0 . 3013E-02  0 . 9671E-03  0.3327E-02 

0 . 73C 

1  c 

0.050 

0.448 

0.9067E+0 4 

0.4472E-02 

0. 1118E+03 
-0.3307E-02 

0. 3180E+04 
-0 . 4118E-03 

0. 7769E+04  0.1410E+04  -0 . 4478E*C4 

0 . 3764E-03  0.7881E-03  -0.7779E-02 

0 . 6869E-03  0.4243E-02  0.3327E-02 

0.82! 

:e 

:  .is: 

0.050 

0 . 1 1 53E+05 
0.4779E-02 

0.2015E+04 

-0.3594E-02 

-0 . 5588E+04 

-0 . 454  9E-04 

0.1236E-05  0 . 11 83E+04  0.4755E-04 
0.6152E-03  0.5697E-03  0.8372E-02 

0. 1836E-02  0 . 2217E-02  0.3824E-02 

C .  78'- 

i*? 

0.151 

0.050 

0.111 1E+05 
0.4779E-02 

0 . 1604E+04 
-0.3594E-02 

0 . 5127E+04 

-0. 4549E-04 

0 . 1149E+05  0 . 12  32E+04  -0.475SE+04 
0.5697E-03  0.6152E-03  -0.8372E-02 

0 . 1574E-02  0 . 2585E-02  0.3824E-02 

0.79! 

16 

c  ,;5; 

0.149 

0 . 12  35E+05 

0 . 4802E-02 

0 . 284  7E+04 

-0. 3558E-02 

-0 . 6414E+04 

-0 . 1564E-03 

0 . 1401E+05  0.1183E+04  0.47S0E-04 
0.7000E-03  0 . 54  36E-03  0.8360E-02 

0.237  6E-02  0.2018E-02  0.3814E-02 

0 . 82C 

^  o 

*  c.  ' 

C.  149 

0 . 1 093E+05 

0 . 4  802E-02 
* 

0 . 1 4  32E+04 

-0 . 3558E-02 

0 . 4  830E*  04 

-0 . 1564E-03 

0 . 1 101E+05  0 . 1352E+04  -0.4750E»04 

0 . 54  36E-03  0.7000E-03  -0.83E0E-C2 

C . I 433E-02  0.3347E-02  0.3814E-C2 

0.85! 
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Appendix  B.  Equipment  List 


Strain  Gages: 

Electrix  Industries:  Lombard  IL 

High  Elongation  Rosettes  (Specially  manufactured) 

Part  Number:  PAHE  -  3  -  125RB  -350  LEN 
Part  Number:  PAHE  -  3  -  062RB  -350  LEN 

Strain  Gage  Adhesive: 

Micro  Measurements:  Raliegh,  NC 
M-Bond  610  High  Elongation 

Recommend  use  as  soon  as  possible  after  application. 

Tab  Adhesive: 

Scotch  3M  Structural  Adhesive  AF  -  163-2 

Miscellaneous  Equipment: 

Instron  20  kip  (thousand  pounds)  universal  test  machine 
Voltmeter 

Summagraphics  Digitizer 

C-Scan  and  Gr/PEEK  Data  sheet  attached.  The  next  four  pages  are  C- 
Scans  of  the  panels  used.  The  last  page  is  the  order  form  describing 
exactly  how  the  material  was  produced. 
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